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ABSTRACT: Starting in 1989, and continuing through the 1990s, high-energy physics witnessed a flowering 
of precision measurements in general and tests of the standard model in particular, led by e'^e~ collider exper- 
iments operating at the resonance. Key contributions to this work came from the SLD collaboration at the 
SLAC Linear Collider. By exploiting the unique capabilities of this pioneering accelerator and the SLD detector, 
including a polarized electron beam, exceptionally small beam dimensions, and a CCD pixel vertex detector, SLD 
produced a broad array of electroweak, heavy-flavor, and QCD measurements. Many of these results are one 
of a kind or represent the world's standard in precision. This article reviews the highlights of the SLD physics 
program, with an eye toward associated advances in experimental technique, and the contribution of these mea- 
surements to our dramatically improved present understanding of the standard model and its possible extensions. 
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During the 1990s, experimental particle physics underwent a quiet revolution. Tests of elec- 
troweak physics in the standard model (^) to unprecedented precision, in some cases at the 
part-per-mil level or better, established the standard model as one of the most successful phys- 
ical theories ever devised. At the center of this activity, starting in 1989, were the experimental 
programs at the Z° boson resonance in e^e colliders at the Stanford Linear Accelerator Center, 
SLAC [using the SLAC Linear Collider (SLC)], and at the European particle physics laboratory, 
CERN [home to the Large Electron Positron (LEP) storage ring]. We review highlights of the 
physics program at the SLC, where a longitudinally polarized electron beam proved to be a 
particularly powerful tool for physics and was critical to an incisive experimental program. 

The weak neutral current, and hence the existence of the boson, is perhaps the most 
significant prediction of the electroweak standard model. Maximal parity violation in charged- 
current weak interactions, implying that the W boson couples only to left-handed fermions, was 
well established and was incorporated by the architects of the standard model. But unlike the 
charged W boson, the neutral Z^ boson couples with both left- and right-handed fermions (with 
different strengths), and the experimental determination of these couplings reveals much about 
the details of the theory. In addition, left-handed and right-handed fermions carry different 
values of a new quantum number, weak isospin, with value Iw = 1/2 for left-handed fermions 
and Iw = for right-handed fermions. The left-handed fermions are grouped into weak-isospin 
doublets, where the upper members (up, charm, and top quarks, and the three neutrinos in 
the lepton case) carry a weak isospin third component {Iw)3 = +1/2, and the lower members 
(down, strange, and bottom quarks, and the charged leptons) carry the value {Iw)3 = —1/2. All 
right-handed fermions are weak-isospin singlets. A consequence of the model is that left- and 
right-handed fermions behave as different particles, underlining the pivotal role of helicity. The 
relative strengths of the left- and right-handed couplings to the Z^ depend on the electroweak 
mixing parameter, sin^ 9w, as follows: 

9L = (%)3 - qsin^ 6w and qr = qsiv? Ow, (1) 

where q is the electric charge. It is useful to define the linear combinations given by 

9v=9L-gR and 9A = gL + gR, (2) 

where gy and qa then appear as the vector and axial- vector couplings, respectively, in the 
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neutral-current Lagrangian density, 



£nc ~ V'(fi'y - QAi^h'^'^Z^. 



(3) 



Here il) and are the fermion and boson fields. Note that the left-handed and right-handed 
terms are included, since {gv - QAl^) = {gv + 9a)\{1 - 7^) + {9V - gA)\{l + 7^)> explicitly 
exhibiting the helicity projection operators. The utility of a longitudinally polarized e~ beam 
for the study of the weak neutral current is clear: With a polarized beam, one can separately 
generate the two different "particles" and e^, and directly measure their couplings to the Z° 
boson. 

The relatively large peak cross section (30.5 nb for hadronic final states) at the pole, 
the clean and well-controlled e"^e~ environment, and the universality of Z^ boson interactions 
(the Z^ couples to all known fermions with comparable strength) combine to make an ideal 
environment for a broad experimental program. As we will demonstrate, when beam polarization 
is available as it was at the SLC, one has a perfect laboratory for testing the standard model. 

At the time of the 1988 International Conference on High Energy Physics direct and 
indirect data led to an uncertainty on the Z^ mass of 0.7 GeV, the limit on the number of 
light Dirac neutrinos was Njj < 6.7 at a 90% confidence level (CL), the error on the electroweak 
mixing parameter of the standard model was ±0.0048 from a combination of all available data, 
the allowed range for the top quark mass was about 50 to 200 GeV, and there was essentially 
no constraint at all on the standard-model Higgs boson mass (~ 5 GeV < niH < ~ 1 TeV). 

Contrast this situation with the present state of affairs (|^): thanks to LEP precision mea- 
surements of the Z^ resonance, the Z^ mass is known to 2.1 MeV and the Z^ "invisible" width 
determines Ni, = 2.9841 it 0.0083; the precise measurement of sin^ 9^, from the SLC Large De- 
tector (SLD) experiment at SLAC and from the LEP experiments at CERN, comes with an error 
of ±0.00017. Electroweak data can now be used to predict the top quark mass with remarkable 
precision (mj = 18lig^ GeV). In addition, the discovery of the top quark at Fermilab in 1994 
and the subsequent determination that rrit = 174.3 ±5.1 GeV (^, in excellent agreement with 
the current standard-model prediction, has greatly improved the predictive power of the theory. 
The standard-model Higgs boson mass is now constrained by the worldwide electroweak data 
to be < 195 GeV at a 95% CL, where Z'^-pole weak mixing angle results are the leading con- 
tributors. Direct searches performed at the higher-energy LEP-II machine constrain this mass 
to be > 114.1 GeV (|5|) — hence, the allowed window for the standard-model Higgs has shrunk 
dramatically. 
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Table 1: SLD data samples 


for the 


various data 


-taking 


; periods 


Year 1992 


1993 


1994-1995 


1996 


1997-1998 


— > hadrons events lOK 


49K 


94K 


52K 


332K 


Polarization 22% 


63% 


77% 


76% 


73% 



The e^e~ Z^-pole programs also contributed to significant advances in heavy-flavor (bottom 
and charm) physics, in particular at the SLC. The unique features of SLC operation led to 
superlative decay-vertex reconstruction in the SLD detector and to highly efficient identification 
of the short-lived bottom and charm hadrons. Thanks to the polarized electron beam, the SLC 
program provides the only direct measurements of parity violation in the bottom and charm 
neutral-current couplings. Heavy- flavor data from LEP, along with current data on the leptonic 
neutral-current couplings (coming from SLC and LEP), can be combined with the SLD results. 
The resulting uncertainties are 2.4% and 1.0% for the vector and axial-vector bottom couplings, 
and 3.7% and 1.1% for the vector and axial-vector charm couplings. The corresponding results 
from the 1988 International Conference were one to two orders of magnitude less precise; for 
example, about 260% and 13% for bottom vector and axial-vector couplings, respectively. 

The most impressive marriage of all of the special features of physics, linear collider op- 
eration, sophisticated detector instrumentation, and beam polarization is evident in the study 
of Bg-B^ mixing at the SLC. These measurements are sensitive to the details of quark mix- 
ing, which are parameterized by the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix (^) in 
the standard model. These data significantly constrain one of the least-known elements of the 
CKM matrix and thereby contribute to our understanding of CP violation. The B'^-B^ mixing 
measurements were also performed at LEP with comparable precision owing to the much larger 
event totals of these experiments. The 1988 International Conference predated these first suc- 
cessful B'^-B^ mixing results, but it is safe to say that the quality of the initial- and final-state 
heavy- flavor identification and the decay vertex resolutions achieved at the Z^, particularly at 
the SLC, were not anticipated. 

The SLD collaboration collected data at the SLC during the period 1992-1998. Table |^ shows 
a breakdown of the data samples, along with the average electron beam polarization at the 
interaction point. Since 1996, data were recorded with the upgraded vertex detector VXD3 (see 
Section |3.2| ). In what follows, we summarize the highlights of the physics program in precision 
electroweak measurements and heavy- flavor physics. The reader is referred to published work 



6 

for studies of r (|^ and QCD (^) physics at SLD. 



ROWSON, SU & WILLOCQ 



2 THE SLAC LINEAR COLLIDER 

The SLC is, so far, the only e~^e~ Unear colhder. Consequently, the SLC/SLD program was both 
a particle physics experiment and an accelerator R&D project, and key features of these parallel 
activities were closely intertwined. A linear collider is a "single-pass" device — the e"*" and e~ 
bunches are generated anew for every single collision — and hence the operation of the SLC intro- 
duced technical challenges in the control of bunch trajectory, intensity, energy, and focusing not 
common to the more forgiving storage ring design. Commissioning of the machine began in 1987 
and lasted about two years. The SLC luminosity was initially far below expectations, although 
the brief Mark II physics program (1989-1990) (^) produced some important early results be- 
fore being eclipsed by the much higher- luminosity LEP e'^e~ storage ring. A breakthrough for 
the SLC came in 1992, when a longitudinally polarized (~ 25%) electron beam first became 
available. Because of advances in photocathode technology, the electron beam polarization was 
increased to ~ 75% by 1994. Over the years, the luminosity of the machine improved steadily, 
most dramatically during the final run (1997-1998), when a peak of about 3 x lO^o cm ^s ^ was 
reached and roughly 70% of the entire SLD dataset was accumulated. 

The cornerstone of the SLC program at the pole was the polarized electron beam, which 
proved essential not only for precision electroweak measurements but for most heavy-flavor 
physics as well. The basic operation of the SLC is described below, where we highlight features 
most relevant to the SLD physics, particularly polarized beam production and transport. More 
complete descriptions are available in an earlier review (|To| ) and in a more recent technical 



report (11) 



2.1 SLC Operation with Polarized Electron Beams 

Figure || illustrates the basic layout of the SLC. The machine produced e'^e~ collisions at a 
repetition rate of 120 Hz. Longitudinally polarized electrons were produced by illuminating 
a GaAs photocathode with a circularly polarized Ti-Sapphire laser (|l^). Since the advent of 
high-polarization "strained lattice" GaAs photocathodes in 1994 (|T3D, where mechanical strain 
induced in a thin (0.1 fim) GaAs layer lifts an angular- momentum degeneracy in the valence 
band of the material, the average electron polarization at the e^e~ interaction point (slightly 
lower than the value produced at the source) was in the range 73% to 77%. Corresponding values 
were about 22% in 1992 using an unstrained "bulk" GaAs cathode, and 63% in 1993 using a 
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Figure 1: Layout of the SLC 



thicker (0.3 fim) strained layer cathode design. The electron helicity was chosen randomly pulse 
to pulse at the machine repetition rate of 120 Hz by controlling the circular polarization of 
the source laser. This pulse randomization minimized unwanted correlations between the pulse 
helicity and other periodic phenomena in the accelerator, as well as in the SLD and polarimeter 
data. On each SLC cycle, two closely spaced (61 ns) electron bunches were produced, the second 
of which was used to produce the positrons. 

The electron spin orientation was longitudinal at the source and remained longitudinal until 
the transport to the damping ring (DR). In the linac-to-ring (LTR) transport line, the electron 
spins precessed in the dipole magnets, where the spin precession angle is given in terms of the 
anomalous magnetic moment g: 



By design, the bend angle ^bend resulted in transverse spin orientation at the entrance to the LTR 
spin rotator magnet. This superconducting solenoid magnet was used to rotate the polarization 
about the beam direction into the vertical orientation for storage in the DR. This was necessary 
because any horizontal spin components precessed rapidly and net horizontal polarization com- 
pletely dissipated during the 8.3 ms (1/120 s) storage time owing to energy spread in the bunch. 
In addition, the polarized electron bunches could be stored in one of two possible configurations 
by reversing the LTR spin rotator solenoid magnet. These reversals, typically done every three 
months, were useful for identifying and minimizing the small [0(10^^)] polarization asymmetries 
produced at the source. 

The product of the horizontal and vertical emittances of the electron (positron) bunches, a 
measure of the transverse phase-space of the beam, was reduced in the DR by a factor of roughly 




(4) 
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100 (1000). One of two stored positron bunches was extracted from the positron DR and led 
the two electron bunches down the linac, the bunches spaced by 61 ns. The second electron 
"scavenger" bunch was diverted to the positron target after 30 GeV of acceleration. Positrons 
produced at the target were accelerated to 200 MeV for transport via the positron return line 
to the positron DR, where they were stored for 16.6 ms. The remaining e~ and e"*" bunches 
continued down the linac and reached ~ 46.5 GeV. They were then transported through the arcs 
(energy loss of about 1 GeV occured in the arcs owing to synchrotron radiation) and final focus 
to the interaction point (IP). Following the IP, the bunches entered the extraction lines where 
their energies were measured by precision spectrometers, after which they were transported to 
beam dumps. 

After leaving the DR, the electrons were accelerated in the linac with their polarization ori- 
ented vertically, and brought into the SLC north arc. The electron spin orientation was manip- 
ulated in the north arc (see Figure |l]) using controllable betatron oscillations (known as "spin 
bumps") to achieve longitudinal polarization at the IP. This was possible because the betatron 
phase advance closely matched the spin precession (1080° and 1085°, respectively) in each of the 
23 bending-magnet assemblies ( "acromats" ) used in the arc — and hence the north arc operated 
close to a spin-tune resonance. As a result, excepting 1992 running, the two additional SLC spin 
rotator solenoids were not necessary for spin orientation and were used only in a series of spe- 
cialized polarization experiments. This unanticipated simplification had an additional benefit. 
In order to achieve higher luminosity, starting in 1993 the SLC was operated with "flat beams" 
{ax > (Jy), a setup that would have been precluded by the spin rotators because x-y coupling in 
the ring-to-linac spin rotator would have spatially rotated the electron bunches. 



2.2 Energy Spectrometers 



The SLC employed a pair of energy spectrometers (14) located in the electron and positron 
extraction lines (see Figure |^). Beam deflection, and therefore beam energy, was measured 
from the spatial separation of synchrotron radiation emitted by the beam in deflector magnets 
located before and after a precision dipole magnet. These two spectrometers together nominally 
measured the e~^e~ collision energy with a precision of 20 MeV. The importance of the energy 



measurement and some details concerning spectrometer precision are discussed in Section 5.2. 
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Figure 2: (a) Horizontal {(ix) and vertical {ay) beam spot sizes, as well as the cross-sectional 
area {ax ■ ay) of the beam as a function of time. The original design value for the area is also 
shown. (6) Luminosity enhancement due to beam-beam disruption. 

2.3 Luminosity Improvements 

The luminosity C of a linear collider is given by 



J- - -. ^d, 

ATraxay 



(5) 



where A^^ are the number of positrons and electrons in each bunch at the IP, / is the collision 
repetition rate, and ax^y are the average horizontal and vertical beam sizes. The additional 
parameter is a disruption enhancement factor that is a function of and the transverse 
and longitudinal beam dimensions, and for moderate beam intensities is equal to one. 

The route to high luminosity at the SLC was somewhat different than originally conceived, 
with a greater emphasis on reduced beam size (design performance was exceeded by about a 
factor of three) because the maximum bunch size was limited to about 4 x 10^'^ due to wakefield 
effects and DR instabilities (7 x 10^'^ was indicated in the initial plans). In the final focus, 
the e~^e~ transverse bunch dimensions were demagnified to about 1.5 in x and 0.7 /um in 
y for collision at the IP. Figure |^ shows the historical trend in the size of the SLC beams. 
A consequence of the very small SLC beam sizes at the IP was that a significant luminosity 
enhancement due to beam-beam disruption {H^ > 1) was observed at the highest luminosities. 
From a comparison between the luminosity measured with the SLD luminosity monitor and that 
computed from Equation |^ with Hd = 1, enhancements of up to a factor of two were seen, as 
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shown in Figure |2^. 

Though not directly relevant to our discussion of the SLD physics program, one important fea- 
ture of SLC operation deserves mention. Producing and maintaining the small beam emittances 
required for these very small spot sizes, and sustaining their collisions, all on a pulse-to-pulse 
basis, was challenging. An important lesson was that a linear collider is an inherently unstable 
machine and extensive feedback systems are indispensable. In particular, the importance of fast 
feedback (at the 120 Hz machine repetition rate) became apparent. By the end of the SLD pro- 
gram in 1998, the SLC featured over 50 feedback systems controlling 250 machine parameters. 
These features included intensity feedback at the source, energy feedback, linac orbit feedback, 
and an elaborate collision optimization system that modulated five final focus corrections for 
each beam while monitoring a direct measure of luminosity (beamstrahlung) . 

2.4 Impact on the SLD Physics Program 

Three features of the SLC contributed to the unique capabilities of the SLD physics program. 
Most important was the polarized electron beam. The SLC achieved a precise determination 
of the weak mixing angle with more than 30 times the sensitivity per event of any compet- 
ing method at LEP and with smaller systematic error. The beam also enabled a number of 
the world's best measurements in heavy-quark electroweak physics, partly through the greatly 
enhanced statistical power and simplicity afforded by polarization. Second, the small size of 
the luminous region at the collision point, smaller in the transverse {x,y) and longitudinal (z) 
dimensions than at LEP by factors of typically 100, 4, and 10, respectively, allowed for excellent 
localization of the primary vertex in reconstructed decays. Finally, the low repetition rate of 
the SLC (120 Hz compared with 45 kHz at LEP), a liability insofar as luminosity is concerned, 
permitted the use of charge-coupled device (CCD) sensors with inherently long readout times 
in the SLD vertex detector, a device that subsequently set the world standard for precision in 
decay-length reconstruction. All these issues are addressed in more detail in our discussion of 
physics results. 

3 THE SLD DETECTOR 

3. 1 Overall Description 

The SLD was a general-purpose device (^) covering most of 47r sterradians (sr) . It was optimized 
for physics at the Z^ resonance, i.e., it was designed for the reconstruction and identification of 
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Figure 3: Side view of one quadrant of the SLD detector. 



particles in a momentum range of several hundred MeV/c up to nearly 50 GeV/ c. The detector 
consisted of multiple concentric layers (see Figure In the following, we refer to a right- 
handed coordinate system with a 2:-axis along the e"*" beam direction and a y-axis along the 
vertical direction. 

Precision charged-particle tracking was performed with the CCD pixel Vertex Detector (VXD) 
surrounded by the Central Drift Chamber (CDC). Outside the tracking systems, the Cherenkov 
Ring Imaging Detector (CRID) identified the charged hadrons tt, and p, and the Liquid 
Argon Calorimeter (LAC) provided electromagnetic and hadronic calorimetry as well as electron 
identification (and muon identification to a lesser extent). These systems were contained inside 
a large solenoid producing a 0.6 T axial magnetic field. Finally, the Warm Iron Calorimeter 
(WIC) served as the magnet flux return yoke and identified muons. In the following, we describe 
the calorimetry systems and devote separate sections to the tracking, particle identification, and 
polarimetry systems. 
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The main calorimeter was the lead-hquid argon samphng calorimeter (^) , which satisfied the 
design requirements for hermiticity (it covered 98% of 47r sr), good energy and position reso- 
lution, and good uniformity of response to electromagnetic and hadronic showers. To provide 
information on longitudinal shower development, the LAC was divided radially into two sepa- 
rate electromagnetic sections totaling 21 radiation lengths and two hadronic sections totaling 2.8 
interaction lengths. Energy resolution was measured to be a/E = 0.15/ E(GeV) for electro- 
magnetic showers and estimated to he a/E = 0.60/ ^JE{GeV) for hadronic showers. The LAC 
was placed inside the magnet coil to minimize the amount of material in front of the calorimeter, 
thus achieving good resolution down to low energies. Finally, the LAC was used to trigger the 
data acquisition system and, in conjunction with the drift chamber, to identify electrons. 

Outside the magnet coil stood the WIC (0), a coarsely sampled iron-gas calorimeter operated 
in limited streamer mode, which was designed to catch the remaining ~ 15% of hadronic energy 
leaking out from the LAC. Its main use, however, was as a muon tracking and identification 
device. The material in the WIC amounted to four interaction lengths, which, in addition to the 
2.8 (0.7) interaction lengths from the LAC (magnet coil), gave a total of over seven interaction 
lengths for particles traversing the whole detector from the IP. 

A third calorimeter system (p^), made of tungsten and silicon pads, detected electromagnetic 
showers at small angle with respect to the beam line: 28 < < 68 mrad. The purpose of this 
system was to monitor the SLC luminosity. 



3.2 The Tracking and Vertexing System 
3.2.1 The Central Drift Chamber 

The CDC (19) had a cylindrical geometry with a length of 1.8m, extending radially from 0.2m 
to 1.0 m. It consisted of ten superlayers, four coaxial to the beam, interleaved with pairs of 
superlayers having stereo angles of ±41 mrad. Each superlayer consisted of cells 50 mm along 
the radius (r) and ~ 59mm wide in azimuth {(p) at the midpoint. Each of the cells contained 
eight sense wire layers. An average spatial resolution of 82 /im was obtained with a gas of 
CO2 : argon : isobutane in the ratio of 75% : 21% : 4%, respectively. The azimuthal and polar 
angle resolutions on the track direction at the inner radius of the CDC were 0.45 mrad and 
3.7 mrad, respectively, for high-momentum tracks. The momentum resolution of combined 
CDC+VXD tracks in the central region (|cos6'| <0.7) was = 0.010 © 0.0024 pr, where is 
the momentum transverse to the beam axis, measured in GeV/c. 
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3.2.2 The CCD Vertex Detector 

The detection of short-hved particles using vertex detectors has played a vital role in the advance 
of our understanding of flavor physics. In the LEP/SLC era in particular, the usage of vertex 
detectors extended beyond the traditional field of heavy-flavor decay physics and became an 
important tool for electroweak and QCD physics, as well as for new particle searches. The low 
beam-crossing rate of the SLC compared with the circular colliders offered an opportunity for 
using CCDs as the basic vertex detector elements. Some of the main advantages over silicon-strip 
vertex detectors employed at LEP are as follows: 

• The SLD vertex detector was a factor of two closer to the beam line, which was crucial for 
reducing track errors due to extrapolation and multiple scattering. 

• The genuine three-dimensional nature of CCD pixel hits led to fewer tracking hit misas- 
signments. In contrast, double-sided silicon-strip detectors yield uncorrelated (j) and z hits, 
which result in a quadratic increase in hit ambiguity. 

• The spatial resolution provided by the SLD CCD pixel segmentation was better by a factor 
of two. 

• The active silicon layer of the CCDs was typically 15 times thinner than that in silicon-strip 
detectors. This feature largely prevented the drastic z-resolution degradation at shallow 
polar angles that is inevitable for silicon strips. It also allowed for thinner detectors to 
reduce multiple scattering. 

The first version of the SLD CCD vertex detector (VXD2) @) was installed in late 1992. Its 
reliable operation and impressive performance inspired the design of an upgraded vertex detector 
(VXD3) (|2l]), with a more optimized geometry, improved coverage, and reduced material. This 
detector was operational starting in 1996 and hence for the majority of the SLD data sample. 
VXDS consisted of 48 "ladders," each containing two 8x1.6 cm custom-made CCDs. The ladders 
were arranged in three cylindrical layers (see Figure ^) at radii of 2.7, 3.8, and 4.8 cm respectively, 
covering |cos0| <0.85 with three layers and |cos0| <0.90 with two layers. The entire detector 
was supported by a beryllium structure and operated at 190 K using nitrogen gas cooling. The 
CCDs were thinned to 200 /um and supported by a beryllium motherboard so that the material 
per layer for normal incidence was 0.40% of a radiation length. The beryllium beam pipe and 
inner gas shell added 0.52% radiation length of material inside the VXD. 

Each VXD3 CCD contained 800x4000 pixels that were 20)Umx20/im in size, and the active 
silicon depth was also 20 /um. VXD3 was continuously read out at a rate of 5 MHz and the entire 
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Figure 4: Layout of VXDS and mechanical support structure showing (a) cross-section view; 
(6) longitudinal section view. 

detector took 200 ms (25 SLC beam crossings) to read out. The electronic noise per pixel was 55 
electrons, compared with the most probable energy loss of 1200 electrons for minimum ionizing 
particles, typically spread over 4-5 pixels. The cluster centroid positions were calculated using 
the measured ionization charge accumulated in each pixel. The average single-hit efficiency was 
over 99.5% for the data taken since 1997; the remaining inefficiency mainly came from a small 
number of faulty electronic channels during limited time periods. The average efficiency for the 
1996 data sample was ~2% worse because of early electronics problems and radiation damage 
effects, which were cured in 1997-1998 by operating at a lower temperature. 

The VXD3 geometry was designed with careful consideration for track-based alignment, with 
CCD overlaps in z on the same ladder and in ^ between adjacent ladders in the same layer. The 
CCD shape distortions were first corrected based on optical survey data. The CCD positions, 
fully described by three translations and three rotations, were subsequently determined using 
a track-based internal alignment. This alignment also included long-range constraints from 
Z° e+e~ and — > n'^n" events. The average single-hit spatial resolution was 3.5 /xm in the 
r-cf) plane for all polar angles and 3.7/xm (5.7 /xm) in the r-z plane for tracks with \cos9\ < 0.6 
(0.6< |cos^| <0.85). The resolution degradation in z at high \cos9\ was partly due to some 
residual errors in the CCD shape distortion corrections. 
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The central beam pipe and VXD assembly was mechanically supported from the CDC. Its 
mechanical coupling to the rest of the beam line was reduced by a factor of ten with a set of 
bellows separating out the central section. This loose coupling provided a very stable relative 
position of the VXD with respect to the CDC, at the level of a few microns, during long periods 
of operation. Global alignment between VXD and CDC was necessary only after major detector 
accesses. 

3.2.3 Tracking Performance 

The tracking algorithm first reconstructed three-dimensional track segments in the VXD with 
three or more hits and pT > 80 MeV/ c. These VXD track segments were then extrapolated to 
the CDC to search for CDC hits to form tracks. On the second pass, the tracks found with 
the remaining CDC hits not used by the first-pass tracks were extrapolated into the VXD to 
search for hits; at least two VXD hits were required. All pattern-recognized tracks were subject 
to a Kalman fit (p2| ) to determine the track parameters, taking into account multiple scattering 
effects. The SLC beam caused a VXD background hit density of 2-7 hits/ cm^ depending on 
the CCD layer and (j) position. This occupancy turned out to be quite manageable thanks to 
the extremely fine granularity of the CCD pixels. For "prompt" tracks in hadronic events, the 
CDC+VXD full tracking efficiency was 96%. Prompt tracks are defined as tracks that originate 
from the IP or from heavy hadron decays with px > 400 MeV/c and | cos9\ < 0.85, and neither 
decayed nor interacted before exiting the CDC. 

Among the prompt tracks with full CDC-I-VXD information, 97.8% had three or more VXD 
hits, and 99.8% had all VXD hits correctly assigned. The VXD track segments without a CDC 
link were used to improve vertex charge determination. Despite the rather short radial span (~ 2 
cm), the charge was correctly determined for 90% of the VXD track segments up to 1.5GeV/c 
(3.5 GeV/c) without (with) a vertex constraint. 

The track impact parameter resolution can be roughly parameterized as ar(f, = 7.7® ^^.^f^2 g 
and arz = 9.6 ^^-^3/2 g /"™) where p is the track momentum, measured in GeV/c. The con- 
stant terms were determined from a study of two-track miss-distance distributions near the IP 
m fi'^fi events, see Figure |. This resolution performance is a factor of two better than 

the best LEP detectors in both the constant and multiple scattering terms. 
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Figure 5: Two-track miss-distance distributions for IJ-'^IJ-" events: (a) r-cf) miss-distance; 

(6) r-z miss-distance. 

3.2.4 The SLC Interaction Point 

To achieve the best identification of secondary vertices from heavy-flavor decays, the primary 
vertex resolution is crucial. To enable the micron-sized SLC beams to collide, it was essential 
to maintain a stable IP. For this reason, the average beam position in the x-y plane is a very 
good estimate of the event primary vertex. In the e+e^ storage ring experiments, the beam 
spot in y is also fairly narrow (a few to 10 /^m), but the beam size in x is typically 100 /im 
or more, so that event-by-event reconstruction of the primary vertex x position was required. 
This typically resulted in a primary vertex x resolution > 20 ^um for LEP experiments. The 
average SLC beam position in the x-y plane was determined from ~30 consecutive hadronic 
events by fitting all tracks coming from a common origin. The resolution of the beam position 
determination was found to be axy = 3.5 /xm using the impact parameter distribution of tracks 
in Z° fJ-'^fJ-~ events. This resolution was the result of the combined effect of beam jitter and 
tracking resolution. 

The event primary vertex in z was measured event by event owing to the large luminous 
region along the beam axis (~700/im). The z position of the primary vertex was determined 
from the median of the z-coordinate distribution for all tracks consistent with the IP in the x-y 
plane. This simple approach was found to be particularly stable against track impact parameter 
resolution tails. The event primary vertex z resolution for {uds, c, h) events was estimated to be 
(10,11,17) Atm. 
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3.3 The Cherenkov Ring Imaging Detector 

The GRID covered the region |cos0| < 0.68. To perform n/K/p separation over the full mo- 
mentum range of interest (~0.2 to 45 GeV/c), two different radiators were used (p^). The liquid 
CqFi/^ radiator identified particles in the lower portion of the momentum range, whereas the gas 
C5F12 + radiator covered the higher end of the range. Thresholds for emission of Cherenkov 
radiation in the liquid (gas) were 0.17 (2.4) GeV/c for vr^, 0.62 (8.4) GeV/c for K"^, and 1.17 
(16.0) GeV/c for protons. Gherenkov photons emitted in a cone around the charged-particle 
trajectory were focused onto an array of 40 time projection chambers (TPGs) filled with ethane 
and ~ 0.1% tetrakis dimethylamino ethylene (TMAE) which served as the photosensitive agent. 
The three-dimensional location of each photoelectron conversion was recorded in the TPGs and 
the production angle relative to a given charged track was reconstructed. For typical /9 ~ 1 par- 
ticles passing through the liquid (gas) radiator, an average of 12.8 (9.2) photons were detected 
per full ring and the Gherenkov angle resolution was measured to be 16 (4.5) mrad. 

3.4 Polarimetry 

In Gompton scattering of longitudinally polarized electrons from circularly polarized photons, 
the differential cross section in terms of the normalized scattered photon energy fraction x is 
given by 

where is the unpolarized differential cross section, V^y and Ve are the photon and electron 
polarizations, and A(x) is the Gompton asymmetry function. The asymmetry arises from the 
difference between cross sections for parallel and anti-parallel spins > (Jj^xl'i)-, ^-i^d is 

defined in the usual way as the cross-section difference divided by the sum of the two cross 
sections. In a polarimeter, the Gompton-scattered photons or electrons are detected, and the 
requisite instrumental effects are incorporated into a detector response function. The normalized 
convolution of A{x) with ^ and the response function (all functions of x) is the "analyzing 
power." When "P-y and the analyzing power are known, the experimental determination of the 
J = 3/2 to j = 1/2 scattering asymmetry determines Ve and hence the utility of this elementary 
QED process to electron polarimetry. 

The SLD precision Gompton polarimeter detected beam electrons that had been scattered 
by photons from a circularly polarized laser. The choice of a Gompton-scattering polarimeter 
was dictated by the requirements that the device be operated continually while beams were in 
collision and that uncertainties in the physics of the scattering process not be a limiting factor 
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in the systematic error — both are troublesome issues for M0ller scattering instruments because 
of their magnetic foil targets. The pulse-to-pulse controllability of the laser polarization, as well 
as the high polarization (99.9%), are additional advantages over other options. 

Figure ^ illustrates the essential features of the polarimeter: Frequency-doubled (532 nm) 
Nd:YAG laser pulses of 8 ns duration and peak power of typically 25 MW were produced at 
17 Hz, circularly polarized by a linear polarizer and a Pockels cell pair. The laser beam was 
transported to the SLC beam line by four sets of phase-compensating mirror pairs and into 
the vacuum chamber through a reduced-strain quartz window. About 30 m downstream from 
the IP, the laser beam was brought into collision with the outgoing electron beam (10 mrad 
crossing angle) at the Compton interaction point (CIP) and then left the beam pipe through 
a second window to an analysis station. The pair of Pockels cells on the optical bench allowed 
for full control of elliptical polarization and was used to automatically scan the laser beam 
polarization at regular intervals in order to monitor and maximize laser polarization at the CIP. 
Downstream from the CIP, a pair of bend magnets swept out the off-energy Compton-scattered 
electrons (typically of order 1000 per laser pulse) that passed through a thin window out of 
the beam- line vacuum and into a nine-channel (1 cm per channel) transversely segmented gas 
Cherenkov detector. By detecting scattered electrons with a threshold Cherenkov device, the 
large beamstrahlung backgrounds in the SLC environment were dramatically reduced. 

The minimum-energy 17.4 GeV electrons, corresponding to full backscattering, generally fell 
into the seventh detector channel. At this point in the electron spectrum, known as the Comp- 
ton edge, the energy-dependent Compton asymmetry function reached its maximum value of 
0.748. Small deviations from the theoretical asymmetry function (of order 1% in channel 7) 
were modeled using an electromagnetic shower simulation ( p4[ ) in a detailed model of the de- 
tector geometry and the magnetic spectrometer (25). The detector was mounted on a movable 
platform and the Compton edge was scanned across several channels at regular intervals in order 
to monitor its location and to experimentally constrain the detector channel response functions. 
This procedure was essential for the precise determination of the analyzing power of the im- 
portant outer channels. Figure |^ presents data showing the corrected Compton asymmetry, as 
well as the magnitude of the correction, as a function of position and scattered electron energy. 
There is good agreement between the corrected asymmetry and the data in each channel. 

Starting in 1996, two additional polarimeters (26), which detected the Compton-scattered 
photons and which were operated in the absence of positron beam due to their sensitivity to 
beamstrahlung backgrounds, were used to verify the precision polarimeter calibration. These 
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Figure 6: A conceptual diagram of the SLD Compton polarimeter. 



two devices were of different design (one employed a gas Cherenkov detector and the other a 
quartz-fiber calorimeter) with different systematic errors, and they had in common with the 
primary electron polarimeter only the instrumental errors due to the polarized laser. The cross- 
check provided by these photon detectors was used to establish a calibration uncertainty of 0.4%. 
Table | summarizes the systematic uncertainties in the polarization measurement. During the 
period 1992-1998, the total fractional systematic error decreased from 2.7% to its final value of 
0.50%, with the most significant reductions coming from greatly improved understanding of the 
laser polarization and the Cherenkov detector nonlinearity. The dominant error is now due to 
the analyzing power calibration discussed above. 

The polarimeter result was corrected for higher order QED and accelerator-related effects (a 
total of -0.22±0.15% for 1997-1998 data). The higher-order QED offset was very smah and well- 
determined (—0.1%) (|27|) . The primary accelerator-related effect arose from energy-polarization 
correlations that caused the average beam polarization measured by the Compton polarimeter to 
differ slightly from the luminosity-weighted average beam polarization at the IP. In 1994-1998, 
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Figure 7: Compton scattering asymmetry as a function of channel position. The horizontal axis 
gives the distance in mm from the center of the detector channels (1 cm wide each) to the path 
of the hypothetical undeflected electron beam (neutral-beam line). The inset shows the seven 
inner detector channels, sized to match the horizontal scale. The per channel data are plotted as 
open circles, and the corrected asymmetry function is the solid curve. The size of the correction 
from the theoretical QED calculation is indicated by the dashed curve. 



a number of measures in the operation of the SLC and in monitoring procedures (smaller and 
better-determined beam energy spread and polarization energy dependence) reduced the size of 
this chromaticity correction and its associated error from its value of 1.1 it 1.7% when it was 
first observed in 1993 to below 0.2%. An effect of comparable magnitude arose from the small 
precession of the electron spin in the final focusing elements between the SLC IP and the CIP. 
The contribution of depolarization during collision was determined to be negligible, as expected, 
by comparing polarimeter data taken with and without beams in collision. All effects combined 
yielded a correction with the uncertainty given in Table ^. 
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Table 2: Compton polarimeter systematic errors on the beam polarization and, in italics, the 



total accelerator-related systematic error 

Uncertainty (%) SVefPe 

Laser polarization 0.10 

Detector linearity 0.20 

Analyzing power calibration 0.40 

Electronic noise 0.20 

Total polarimeter uncertainty 0.50 
Chromaticity and interaction point corrections 0.15 



Table ^ gives the fully corrected, luminosity-weighted, average polarizations corresponding to 
each of the SLD runs. Improvements in GaAs photocathode performance are evident in the 
1993 run (first use of a strained-lattice material) and the 1994-1995 run (in which the active 
layer was three times thinner) . Changes in the achieved polarization in later years mainly reflect 
variation in photocathode manufacture. 

A number of experiments and redundant systems were used to verify the high-precision po- 
larimeter. Most important were the following: 

• Moderate precision M0ller and Mott polarimeters confirmed the high-precision 
Compton polarimeter result to ~ 3% (1993 1995), and gamma polarimeters 
were operated in parallel with the electron detector polarimeter (1996 1998). 

M0ller polarimeters located at the end of the SLAC linac and in the SLC electron extraction 
line were used to cross-check the Compton polarimeter. The perils of using a less reliable 
method to test a precision device were apparent when large corrections for atomic electron 
momentum effects in the M0ller target were discovered (^) and accounted for, after which 
good agreement was obtained. In addition, a less direct comparison was provided by Mott 



polarimeter bench tests of the GaAs photocathodes (29). As described above, during the 



Table 3: Luminosity-weighted average polarization values for all SLD data 



1992 


1993 


1994-1995 


1996 


1997-1998 


0.224 


0.630 


0.7723 


0.7616 


0.7292 


± 0.006 


±0.011 


± 0.0052 


± 0.0040 


± 0.0038 
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last years of the program, two photon detectors were used to cross-check the primary 
electron polarimeter and to experimentally establish calibration uncertainties. 

• SLC arc spin transport was extensively studied (1993 1998) and was frequently 
monitored and adjusted. A series of experiments studied the beam polarization re- 
ported by the Compton polarimeter as a function of beam energy, beam energy spread, 
and beam trajectory in the SLC arcs (pO|). Two spin rotators (in the linac and in the ring- 
to-linac return line) were scanned in order to determine the IP polarization maximum. An 
important result of these experiments was the discovery that the SLC arcs operate near a 
spin tune resonance, leading to the advent of spin manipulation via "spin bumps" in the 
SLC arcs mentioned above. This procedure eliminated the need for these spin rotators 
and allowed one to minimize the spin chromaticity (dVe/dE), which reduced the resulting 
polarization correction from > 1% in 1993 to < 0.2% by 1995. 



4 THE PHYSICS ANALYSIS TOOLS 



4-1 Monte Carlo Simulation 



The Monte Carlo simulation for — > hadrons events employed the Jet set 7.4 generator (31). 
For the B meson decay simulation, the CLEO QQ generator was adopted, with adjustments 
to improve the agreement with inclusive particle production measurements from ARGUS and 



CLEO. A detailed simulation of the detector was performed with Geant 3.21 (32). Simulated 
decays were overlaid with signals from events taken on random beam crossings in close time- 
proximity to each recorded Z^ and then processed using the standard data reconstruction. The 
Monte Carlo Z^ hadrons events were produced with statistics exceeding data by factors of 4, 
15, and 22 for uds, cc, and bb events, respectively. 

For many of the measurements described in this article, the systematic errors associated 
with tracking efficiency and resolution uncertainties are important and require a consistent 
evaluation. The uncertainty in the tracking efficiency was evaluated from a comparison between 
data and Monte Carlo for the fraction of good CDC tracks extrapolating close to the IP that 
did not have associated VXD hits. These studies indicated that the Monte Carlo simulation 
overestimated tracking efficiency by ~1.5% on average. A procedure for the random removal 
of tracks in bins of pT, (p, and was used to correct the Monte Carlo for this difference. Most 
analyses quoted the effect of the entire correction as a systematic error. The effect of residual 
detector misalignment on the tracking was estimated from the observed shifts in the track impact 
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parameter distributions as a function of (j) and 6 in both r-(j) and r-z planes. The typical impact 
parameter biases observed in the data were ~ 2.5 /um (~ 5//m) in the r-<j) (r-z) plane. A 
correction procedure was applied so that the Monte Carlo tracks match the mean bias values of 
the data in various {(j), 6) regions. This was a more realistic evaluation of alignment bias effects, 
where tracks passing the same detector region were biased in a correlated manner. The full 
effect of this correction was taken as the systematic error. 

4-2 Vertexing Algorithm 

Conventional identification of heavy-flavor decays with vertex detectors typically relies on either 
counting tracks with large impact parameter significance or vertexing kinematically selected 
decay tracks. SLD took advantage of its CCD vertex detector and the small and stable SLC 
IP to pioneer inclusive topological vertexing methods that provided additional information on 
decay kinematics, charge, and flavor. 

The topological vertexing algorithm ZVTOP ( |33| ) aimed to identify all separated vertices and all 
secondary decay tracks in the event. Well- measured tracks with VXD hits were used to search 
for secondary vertices by looking for areas of high track overlap density in three-dimensional 
coordinate space, taking into account the individual track resolution functions. Because of 
the excellent decay length resolution available at SLD, it was possible in a significant fraction 
of decays to distinguish secondary (6-hadron decay vertex) tracks from tertiary (cascade charm 
decay vertex) tracks, instead of lumping all secondary tracks into a single vertex, as was typically 
done by other experiments. The vertexing algorithm was therefore designed to allow the search 
for multiple vertices associated with one 6-hadron decay. 

The vertexing was done in each hemisphere or jet separately. The b hadrons produced in 
decays have an average momentum of ~30GeV/c, corresponding to 7/5 ~6. This large boost 
distributes the tertiary vertices along the 6-hadron flight direction. The original algorithm, as 
described by Jackson (|3^), simply used a 50/um-radius cylindrical volume around the jet or 
hemisphere cLXlS clS cl preferred secondary vertex volume to suppress fake vertices. An improved 
version of ZVTOP, used by some of the analyses in this article, deployed a "ghost" track to 
best represent the 6-hadron flight direction. This was achieved by finding the direction that 
minimized the vertex sum. for the vertices between the ghost track and all other tracks in the 
jet. It should be noted that once the ghost track was removed, one could be left with single-track 
"vertices," which is in fact a fairly common decay topology (e.g., one-prong decays). This 
algorithm found at least one secondary vertex in 73% (29%) of the hemispheres in bb (cc) events. 
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Figure 8: Variables used in ZVTOP track-vertex assignment. 

Among the b hemispheres that had at least one secondary vertex, two or more secondary vertices 
were found in 30% of them. 

To obtain the best description of a 6 hadron, it is crucial to identify all its decay products. 
Since the basic decay vertex topology was established with well-measured tracks only, an attempt 
was made to incorporate lower-quality tracks or even unlinked VXD track segments. To do so, 
a vertex axis was formed by joining the IP to a "seed" vertex combining both secondary and 
tertiary tracks, and the distance D between the IP and this vertex was computed. For each track, 
the three-dimensional distance of closest approach T to the vertex axis and the longitudinal 
distance L between the IP and the point of closest approach on the vertex axis were calculated 
(see Figure ^). Typical requirements for attaching a track to the seed vertex were T < 1 mm, 
L > 250 /xm, and L/D > 0.25. This track attachment procedure was further improved by means 
of a neural-network algorithm combining the information on T, L, L/D, and track angle with 
respect to the vertex axis. There were on average 3.9 quality secondary tracks in the seed vertex, 
and 0.9 additional tracks were attached per b hemisphere. The VXD track segments were used 
only for vertex charge determination, with an average of 0.2 such segments attached per 6-hadron 
decay. The sum of the identified secondary tracks and VXD track segments corresponds to an 
average of 82% of all prompt 5-hadron decay tracks, with a track assignment purity of 96.8%. 

4-3 b and c Tagging 

Prior to the advent of vertex detectors, events with b or c quarks were typically identified by the 
presence of high-p^^ leptons or fully reconstructed D* decays. However, these tagging methods 
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suffered from relatively low branching ratios as well as limited purity. The development of 
vertex detectors opened a new avenue in heavy-flavor tagging by exploiting the surprisingly 
long 6-hadron lifetime. Limited success in 6-tagging had already been achieved with vertex drift 
chambers at PEP/PETRA (|3^), but 6-tagging developed into a powerful tool in the LEP/SLC 
era only with the advent of silicon vertex detectors. SLD contributed pioneering vertexing 
techniques to fully explore the high resolution offered by the CCD pixel vertex detector and 
elevated inclusive b- and c-tagging to a new level. 

There are various ways of using the basic tagging tools (tracks with large impact parame- 
ters and secondary vertices), which are all quite efficient for 6-tagging. However, it is difficult 
to achieve 5-tag purities > 95% without drastic loss of efficiency because the main remaining 
background consists of irreducible cc events with genuine long decay length. An important 
improvement in 6-tagging was the utilization of the secondary vertex mass as a powerful dis- 
criminator to separate b and c quarks. 

The SLD implementation of the vertex mass tag included additional kinematical information 
to improve the efficiency while still maintaining high purity. The tracks identified as secondary 
by the procedure described in the previous section were used to calculate the raw vertex invariant 
mass (Mraw) for each hemisphere. For vertices with Mraw close to the charm cutoff of 2 GeV, 
the vertex momentum is typically collinear with the vertex axis for vertices in cc events whereas 
vertices in bb events at this low mass typically have large missing transverse momentum (pt) with 
respect to the vertex axis because of missing neutral particles. This information was utilized by 
forming a "pj-corrected mass": 



This minimum-pf correction can be simply derived by boosting the system of identified secondary 
tracks along the vertex axis until their longitudinal momentum sum is zero, and assuming that 
a massless neutral system is recoiling with the missing pt in that frame. To prevent udsc events 
from gaining an artificially large pt due to detector resolution effects, the pt was actually the 
minimum pt allowed by the uncertainties in the primary and secondary vertex positions. Afcorr 
was further restricted to be < 2 x M^aw to reduce the contamination from fake vertices in 
light-quark events. The effectiveness of Afcorr as a 6-tagging variable can be clearly seen in 
Figure ^a. Many analyses using only a basic secondary track classification and requiring a mass 
tag of Mcorr > 2 GeV benefited from a hemisphere 6-tagging efficiency of 53% and a b purity 
of 98%. In this simple scheme, a cut on the raw mass of Mraw > 2 GeV would give a similar 
purity, but the efficiency would be only 35%. The partial correction for missing neutral particles 




(7) 
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Figure 9: Distributions of (a) the pt-corrected mass Mcorr and (6) the neural-net c — b separation 
variable Scb comparing data and Monte Carlo. 

depends critically on the ability to precisely measure the primary and secondary vertex positions 
to ensure that the vertex flight direction was well determined. The fact that so far only SLD 
has effectively applied this underscores the crucial role of the small and stable SLC IP and the 
high-precision vertexing achieved with the CCD vertex detector. 

Although the requirement Mcorr > 2 GeV cleanly selected the high-mass b hemispheres, further 
kinematical information can be used to achieve good separation between b and c in the mass 
region below 2 GeV. A neural-network algorithm was constructed to process information on the 
secondary vertex momentum, track multiplicity, and decay length, in addition to the p^-corrected 
mass, leading to a single c — h separation variable S^t as shown in Figure |9^. The b tags were 
given by large values of Scb, whereas a clean c tag also emerged at small Set- A typical b tag 
requiring Set > 0.75 gave a hemisphere 6-tagging efficiency of 65=62% and purity of nf,=98.3%. 

We not only enhanced the 6-tag performance but also simultaneously obtained a c tag with 
a typical cut of Set < 0.30, giving ec=18% and nc=84%. Given the shorter lifetime and much 
less distinctive kinematical properties of cc events, this is also a very impressive performance, 
which is approaching the earlier era of 6-tagging performance. It should be noted that the uds 
events amounted to only a very small fraction of the events containing a secondary vertex, and 
they were concentrated near Scb = 0.5. Furthermore, a large fraction of the uds background was 
actually due to g ^ bb and g ^ cc production in uds events. The success of this inclusive charm 
tag is unique to SLD. 
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4.4 Hadron Identification 

SLD was one of two resonance experiments with a Cherenkov ring imaging system for hadron 
identification. The primary goal of the GRID was to identify charged pions, kaons, and protons, 
but the system proved valuable in lepton identification also (see next section). Hadron identifica- 
tion at SLD relied exclusively on the GRID, i.e., no dE/dx information from the CDG was used. 
Tracks were extrapolated outward from the end of the drift chamber through the liquid and gas 
radiators. The TPC hit positions were then transformed into Gherenkov angles with respect to 
the track trajectory for each of the radiators. Finally, a global likelihood approach was used to 
identify the charged particles. A likelihood Ci was computed for each particle hypothesis (i = 
e, n, IT, K, or p). The likelihood made use of the expected number of photons and Cherenkov 
angle, which both depend on the chargcd-particle momentum and the mass hypothesis chosen, 
the number and angles of reconstructed photons, and a background term. The latter incorpo- 
rated information about background from random hits as well as the contribution from genuine 
Gherenkov photons originating from other charged particles in the event. The identification pro- 
ceeded by cutting on differences between the logarithms of these likelihoods (the log-likelihoods 
from the two radiators were simply added to form an overall likelihood). 

Kaon identification is very important in heavy-flavor physics, and it is preferable to keep as 
high an efficiency as possible by focusing mostly on pion rejection. Proton rejection was difficult 
in the momentum range 3 < p < 9 GeV/ c because both kaons and protons were below Gherenkov 
threshold in the gas radiator and were thus primarily identified by the absence of Cherenkov 
radiation. Pion rejection was achieved by requiring log Ck — log/^Tr > 3 (> 5) for p < 2.5 
(> 2.5) GeV/c, corresponding to the momentum region below (above) threshold for charged pion 
radiation in the gas radiator. Loose proton rejection was obtained with log — log^Cp > — 1. 
Other cuts were also applied to guarantee that the detector was fully operational and that the 
quality of the track extrapolation was high. For kaons with p > 0.8 GeV/c and |cos^| < 0.68, 
the selection had an efficiency between 40 and 55%, depending on the momentum of the track. 
Using clean pion samples from Kg tt+tt" decays, the tt ^ K misidentification rate was 
measured to be nearly constant at 2.5% for 0.8 < p < 2.5 GeV/c, rising to about 10% between 

2.5 and 4 GeV/c and remaining ~10% for p > 4 GeV/c. 

4.5 Lepton Identification 

Like hadron identification, lepton identification required that charged tracks be extrapolated 
outward from the drift chamber to the LAG and WIG subsystems. The extrapolation transported 
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the track parameters as well as the error matrix, taking multiple scattering into account. 

Electron identification was performed in the range | cos9\ < 0.72 for tracks with p > 1 GeV/c 
and relied primarily on information from the LAC. A special feature of the identification was 
the use of the GRID to remove charged hadron contamination. A neural network, tuned to the 
Monte Carlo simulation and cross-checked with data, optimized the electron selection based on 
the LAC and GRID variables. The selection efficiency was estimated using the Monte Carlo 
simulation to be 64%, whereas the purity was 64% for p > 2 GeV/c. A pion misidentification 
rate of approximately L0% was measured with a clean sample of pion tracks from Kg -jt'^'k~ 
decays, in good agreement with the simulation. Electrons from photon conversion were identified 
and rejected with a 73% efficiency. 

Muon identification was performed in the range | cos^| < 0.70. Extrapolated charged tracks 
with p > 2 GeV/c were matched with hit patterns in the WIG. Guts on the quality of the hit 
pattern and CDG/WIC match, as well as the penetration depth, provided a sample of muons 
with 81% efficiency and 68% purity, as determined in the simulation. This muon selection 
also took advantage of the GRID to reject approximately half of the charged kaon and proton 
contamination in most of the momentum range, and roughly a third of all charged pions with 
p < 6 GeV/c (with an efficiency loss of only 5%). Information from the pattern of energy 
deposition in the LAG was also used in the identification and was especially useful at larger 
|cos^|. Pion misidentification rates of about 0.3% were estimated from — tt+tt" and t"^ 
decays in the data, in reasonable agreement with the simulation. 

5 ELECTROWEAK PHYSICS: BOSON COUPLINGS 

5. 1 Introduction 

In Section 1, we mentioned that in the standard model the Z f f couplings depend on the weak- 
isospin of the fermions and on a single parameter, sin^ 6y/, but we did not discuss the origin of 
this parameter. Recall that in the standard model, the weak isotriplet and the isosinglet 
gauge fields, with gauge couplings g and g' , arc mixed by electroweak symmetry breaking due to 
the finite vacuum expectation value of the Higgs scalar field. The charged fields and and 
neutral fields and combine linearly into the physical charged and W~ and the neutral 
photon and Z^ gauge bosons. This change of basis is parameterized by a weak mixing angle 
Ow, given by g' = gtai,n{6w), and the masses of the physical gauge bosons (and the fermions) 
emerge as by-products of electroweak symmetry breaking (the "Higgs mechanism" ) . 
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Electroweak tests of the standard model reached an important turning point once the 
boson mass was determined at LEP to a precision of two parts in 10^. The measurement of Mz 
provides a third precision constant, which together with the Fermi constant Gp (constrained by 
muon decay) and the fine structure constant a (evaluated at = M^) is sufficient to determine 
the three universal parameters of the electroweak standard model: the SU{2)i^ x U (1) couplings 
g and g' , and the vacuum expectation value of the Higgs field. The couplings of fermions to the 
boson are, by virtue of weak mixing, a function of 9w — hence, their determination provides 
a fundamental test of electroweak symmetry breaking and, if sufficiently precise, of higher-order 
corrections. The electroweak measurements made by the SLD collaboration can generally be 
described as measurements of the fermion-to-.Z^'^ couplings. 

The differential cross section for e^e'^ —fZ^—f ff is expressed as 

= (1 - VeAe)il + COS^ 9) + 2 COS e{Ae-Ve)Af, (8) 

a cos t> 

where cos is the cosine of the angle between the final-state fermion / and the incident electron 
directions, Ve is the electron beam longitudinal polarization, and and Af are the asymmetry 
parameters for the initial- and final-state fermions, respectively. The asymmetry parameter for 
a given fermion represents the extent of parity violation at the Z'' — > // vertex and is defined 
as 

_ 2gvgA _ gj-gji 

The parameter Af can be isolated by the measurement of various cross-section asymmetries, 
which is experimentally attractive because systematic effects are minimized. In addition, gy+g\ 
(or equivalently, g£ -|- g]^) is determined from the Z^-decay partial widths, which are normalized 
by the hadronic decay partial width to control systematic effects. 

The simplest of the asymmetries is the left-right cross-section asymmetry 

^L = ^^=^e, (10) 

for which all angular dependence and all dependence on the final state cancej^. As a result, 

Alr is a particularly robust quantity, with smaller systematic effects than any other asymmetry. 

This asymmetry provides a direct measurement of the coupling between the Z^ and the e~^e~ 

initial state, and provides, as we shall see, the best sensitivity to sin^ dw- 

^The dependence of ^lr on the final state couplings and polar (and azimuthal) angle completely vanishes, 
provided that the efficiency for detecting a fermion at some polar angle (with respect to the electron direction) 
is equal to the efficiency for detecting an anti-fermion at the same polar angle. This condition is satisfied by the 
SLD detector. 
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Table 4: The approximate magnitude of the various fermion-to-Z'^ couphng parameters, for 
sin2 = 0.23 



e,n,T -0.27 -0.23 0.05 0.15 -7.9 

u,c 0.35 0.15 0.17 0.67 -3.5 

±s,b -0.42 -0.08 0.22 0.94 -0.6 



Asymmetries that retain angular information are sensitive to the final-state couplings. For 
~^ ff decays, the forward-backward asymmetry, measured at LEP for lepton and heavy 
quark final states, can be expressed in terms of z = cos 9 as 

aHz)-af{-z) _ 2z 

The asymmetry for fermions is a composite observable sensitive to both the initial-state 
and the final-state Af. For example, from the magnitudes of the parity violation parameters Af 
as listed in Table ^, it can be seen in the case of the b quark that the large value of Ai^ makes A^^ 
particularly sensitive to A^ (and hence to sin'^^iy). In general, the final-state Zff asymmetry 
parameter Af can be deduced by taking the lepton coupling ^e's from other measurements such 
as the r polarization, the lepton pair forward-backward asymmetries, and j4lr. Note that Ap^ 
is subject to systematic errors in the determination of detector acceptance and efficiency. 

At the SLC, the polarized forward-backward asymmetry can be measured. Due to the factor 
{Af. — Ve) in Equation ^, manipulation of the helicity of the e~ beam (Pg < for left-handed 
electrons) distinguishes two different forward-backward asymmetries. In particular, for a highly 
polarized beam, the forward-backward asymmetry is not only of the opposite sign for left- 
handed beam, the magnitude is also expected to be larger compared with the asymmetry for 
right-handed beam. Separate measurements of the left and right cross sections can be combined 
into the left-right forward-backward asymmetry 

AU{z) = ^-^f-^ ^ ^-^^ ^ ^ = \Ve\Af ^. (12) 

™ k{(z)+a{(-z)] + [4(.)+4(-z)] -^1 + ^^ ^ ^ 

The use of A^-^ eliminates the dependence on A^. and measures Af directly. Compared with 
Apg, benefits from a large gain in statistical power for the determination oi Af. Given an 
SLC electron beam polarization of ~75%, this improvement factor is {Ve/Af.)"^ ~ 25. The effects 
of nonuniformity of the detector acceptance and efficiency cancel to first order for this double 
asymmetry. The SLD collaboration has measured A^-q for bottom, charm, and strange quarks as 
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well as the charged leptons, providing direct measurements of the associated fermion asymmetry 
parameters. With the exceptions of and Ar, these direct measurements are unique to SLD. 

Measurements of partial Z^-decay widths provide information complementing that obtained 
from the asymmetries. Because 

'^^ = r(ZO^hadrons) "^^^ + ^^' ^'^^ 
Rf measures the — > // coupling strength, whereas Af measures the extent of parity violation 
at the —> ff vertex. From another perspective, consider for example the sensitivity of Rf, 
and Ah to the left- and right-handed Zbb couplings: 

6Rb/Rb ~ -3.57(55^-0.65(55^ and 

(5A/A ~ -0.31(55i + 1.72(5(7^. (14) 

We see that Rj, is more sensitive to the possible deviations in the left-handed Zbb coupling, 
and Ab is more sensitive to deviations in the right-handed coupling. The SLD collaboration has 
measured both Rb and Re- 

We have discussed how the well-determined constants Mz, Gp, and a(M^) constrain the 
standard model. However, these data are a sufficient set only at tree level, and to test the 
standard model, including loop effects, at least one additional precise measurement is required. 
In the standard model, aside from the effects of real and virtual photons, radiative corrections 
to the fermion-to-gauge boson couplings are dominated by corrections to the boson propagators. 
These effects are known as oblique corrections because the virtual loop is not directly coupled 



to the initial- or final-state fermions (see Figure 10). Loop corrections to the boson propagators 
are dominated by those due to the heaviest fermion, the top quark, and are proportional to 
(mf — m^)/M^ (note that this factor depends on weak-isospin symmetry breaking in the top- 
bottom doublet). This arises because of the proportionally large Yukawa couplings of the heaviest 
fermions to the Higgs boson. The sensitivity of the radiative corrections to rrit means that 
precise electroweak measurements can be used to constrain nit. Now that nit has been measured 
directly, precision measurements have become sensitive to much smaller effects, such as those 
due to the Higgs boson, which are only logarithmic in the Higgs mass, and those due to new 
physics processes. 



In the standard model, direct vertex corrections (see Figure 11) are small, but they are ex- 
pected to be largest for the heaviest fermions. In particular, the Zbb vertex corrections, which 
act only on the left-handed coupling, are particularly large owing to the large top quark mass 
and the fact that the relevant quark-mixing matrix element {\Vtb\ ~ 1) is large. In general, 
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Figure 10: Feynman diagrams showing radiative corrections affecting the propagator, includ- 
ing contributions from fermions, gauge bosons and the Higgs boson. 

heavy quark partial widths and asymmetry parameters are most sensitive to vertex corrections 
but with different sensitivity to left- and right-handed coupling constants. In contrast, oblique 
corrections are best isolated by measurements of the lepton asymmetries, and this distinction has 
an important consequence. In contrast to the sin^ 9w sensitive lepton asymmetries, the vertex- 
sensitive observables are largely independent of the Higgs boson and top quark masses, so their 
standard-model predictions are relatively unambiguous. Hence, the SLD precision electroweak 
results constitute a complementary set of neutral-current coupling measurements. 



5.2 Lepton Couplings 
5.2.1 The Alr Measurement 



The measurement of the left-right cross-section asymmetry Alr (35,^) at the SLC provides a 
determination of the asymmetry parameter A^. and is presently the most precise single measure- 
ment, with the smallest systematic error, of this quantity. In addition, Ae is particularly sensitive 




f e- 



f e- 



Figure 11: Feynman diagrams showing radiative corrections affecting the coupling to the 
final state fermions. 
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to the effective weak mixing angle witli 6A(, ~ 7.85((5 sin^ (for tlie Apg measurements, tlie 
analogous factor is 1.77). Hence, the most precise available determination of the effective weak 
mixing angle derives from the ^lr measurement. 

In principle, the analysis is straightforward. One counts the numbers of bosons produced by 
left- and right-handed longitudinally polarized electrons (A^l and N^), forms an asymmetry, and 
then divides by the luminosity- weighted e~ beam polarization (the e"*" beam is not polarized): 



The method requires no detailed final-state event identification (e^e~ final-state events are 
removed because of their nonresonant t-channel contributions, as are all other backgrounds not 
due to decay) and is insensitive to all acceptance and efficiency effects. In order to convert 
the j4lr into a determination of the effective weak mixing angle, the result is converted into a 
"Z'^-pole" value by the application of a small (typically about 2.0%) correction for initial-state 
radiation and 7 — Z*^ interference ( |37D : 



This calculation requires accurate and precise knowledge of the luminosity-weighted average 
center-of-mass collision energy Ecm- 

For the most recent data (1997-1998), the small total systematic error of a relative 0.65% is 
dominated by the 0.50% relative systematic error in the determination of the e~ polarization, 
with the second largest error (0.39%) arising from uncertainties in the determination of the 
luminosity-weighted average center-of-mass energy. Some much smaller contributions to the 
systematic error are discussed below. The relative statistical error on Alr from all data is about 



Below, we describe some details of the ^lr measurement and provide some historical context 
for the Alr program at SLC/SLD (1992-1998). 

Systematic Effects and Their Control The Alr measurement is remarkably resistant to 
detector-dependent systematic effects and Monte Carlo modeling uncertainties that are signifi- 
cant issues for most other electroweak precision measurements. By far the dominant systematic 
effects arise from polarimetry and from the determination of the collision energy, rather than 
from any details of the analysis or the operation of SLD. The simple expression given in Equa- 
tion |l^ applies to the ideal case in the absence of additional systematic effects, and as such it is 
a good approximation to better than 0.2%. Nevertheless, systematic left-right asymmetries in 
luminosity, polarization, beam energy, and acceptance, as well as background and positron polar- 



^LR 



A^L-A^R 1 

A^L + iVR {Ve) ■ 



(15) 




(16) 



1.3% 
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ization effects, can be incorporated into an extended expression for the cross-section asymmetry. 
One finds that the measured asymmetry Am is related to by the fohowing expression, which 
incorporates a number of small correction terms in lowest-order approximation: 



{Ve) {Ve) 



/bkg(^m - ^bkg) - Ac+ A^A-p 

(17) 



E,^^^j^AE-A + {Ve)V, 



where {Ve) is the mean luminosity-weighted polarization; /b^g is the background fraction; cr{E) 

is the unpolarized boson cross section at energy E; cr'{E) is the derivative of the cross section 

with respect to E; ^bkg; Ac, A-p, Ae, and A^ are the left-right asymmetrie^ of the residual 

background, the integrated luminosity, the beam polarization, the center-of-mass energy, and the 

product of detector acceptance and efficiency, respectively; and Vp is any longitudinal positron 

polarization, which is assumed to have constant helicity. Because the colliding electron and 

positron bunches were produced on different machine cycles and because the electron helicity 

of each cycle was chosen randomly, any positron helicity arising from the polarization of the 

production electrons was uncorrelated with the electron helicity at the IP. The net effect of 

positron polarization from this process vanished rigorously. However, positron polarization of 

constant helicity would affect the measurement. 

The close ties between the measurement and the SLC accelerator complex are highlighted 

by numerous accelerator-based experiments dedicated to the SLD physics program. Particularly 

important among these experiments was a Z'^-peak scan to calibrate the sum of the beam 

energies. The SLC energy spectrometers were briefly described in Section These devices were 

first operated in their final configuration in 1989 by the Mark II experiment, and the calibration 

of the two precision spectrometer magnets was performed in 1988 (|l^). Their expected precision 

was about ±20 MeV on the measured center-of-mass collision energy Ecm- The importance of 

these devices to the measurement is quantified by the approximate rule of thumb that an 

80 MeV uncertainty in Ecm corresponds to a 1% error on the Z'^-pole asymmetry A^j^. For 

this reason, a Z'^-peak scan was performed in 1998 to calibrate the spectrometers to the LEP 

measurement of the mass. The scan used two optimized off-peak points at -1-0.88 and —0.93 

GeV and a luminosity approximately equivalent to 9000 produced Z^ bosons to reach a statistical 

precision on the peak position of 20 MeV. The results of a complete analysis of systematic effects 

determined an offset of —46 MeV and a total Ecm uncertainty of 29 MeV [the 0.39% uncertainty 
■^The left-right asymmetry for a quantity Q is defined as Aq = {Ql — Qr)/{Ql + Qr) where the subscripts L 
and R refer to the left- and right-handed beams, respectively. 
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on vl^R mentioned above and given in Table |5| (^)] . 

Other examples of accelerator-based experiments include the following: 

• The e bunch helicity transmission was verified by setting up a large current- 
helicity correlation in the SLC, allowing for the use of the LAC to verify data 
synchronization (1992—1993). Although the electron bunch polarization state was 
transmitted via reliable and redundant paths to the SLD detector/polarimeter complex, 
the SLD electroweak group proposed a series of independent tests of the synchronization 
of this information and the SLD event data. In one such test, the laser optics at the SLC 
polarized source were temporarily modified by the addition of a polarizer and quarter-wave 
plate so that photocathode illumination was nulled for one of the two circular polarization 
states. The positron beam was turned off, and the electron beam was delivered to the IP. 
Beam-related background in the LAC was detected, but only for the non-extinct pulses. 
By this means, the expected correlation between helicity and the presence of beam, and 
hence the LAC data stream, was verified (^). 

• Positron polarization was experimentally constrained. In 1998, a dedicated exper- 
iment was performed in order to directly test the expectation that accidental polarization 
of the positron beam was negligible. The positron beam was delivered to the SLAC End 
Station A (ESA), where a M0ller polarimeter was used. Experimental control was assured 
by first delivering the polarized electron beam, and then an unpolarized electron beam 
(sourced from SLAC's thermionic electron gun), to the ESA, confirming polarimeter op- 
eration. In addition, the spin rotator magnet located in the linac was reversed halfway 
through the positron beam running, thereby reversing the sense of polarization at the 
M0ller target and reducing systematic error. The final result verified that e"*" polarization 
was consistent with zero (-0.02 ± 0.07%) (|0|). 

A simple calorimetric event selection in the LAC, supplemented by track multiplicity and 

topology requirements in the CDC, was used to select hadronic decays. For each event 

candidate, energy clusters were reconstructed in the LAC. Selected events were required to 

contain at least 22 GeV of energy observed in the clusters and to have a normalized energy 

imbalance]^ of less than 0.6. The left-right asymmetry associated with final-state e~^e~ events 

is expected to be diluted by the t-channel photon exchange subprocess. Therefore, we excluded 

e~^e~ final states by requiring that each event candidate contain at least four selected CDC 
^The energy imbalance is defined as a normalized vector sum of tlie energy clusters as follows, -Eimb = 

cluster • 
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tracks, with at least two tracks in each hemisphere (defined with respect to the beam axis), or at 
least four tracks in either hemisphere. This track topology requirement excluded Bhabha events 
that contained a reconstructed gamma conversion. Small backgrounds in the Alr data sample 
were due to residual e~^e~ final-state events, and to two-photon interactions, beam-related noise, 
and cosmic rays. For the most recent data (1996-1998), the total background contamination 
was estimated to be < 0.05% for a selection efficiency of (91 ± 1)%. 

The asymmetries in luminosity, polarization, and beam energy were all continually monitored 
using small-angle Bhabha counters, the Compton polarimeter, and energy spectrometer data 
available at the SLC repetition rate of 120 Hz. These asymmetries were limited to approximately 
10~^, 10~^, and 10^^, respectively. The long-term average values of all asymmetries of this 
type were reduced by the roughly trimonthly reversal of the transverse polarization sense in 
the electron damping ring (DR) referred to in Section |2|. The dominant cause of the observed 
asymmetries was the small current asymmetry produced at the SLC polarized source. This 
effect arose because of the source photocathode sensitivity to linear polarized light, together 
with residual linear polarization in the source laser light that was correlated with the light 
helicity. The source current asymmetry was minimized by a polarization control and intensity 
feedback system (starting in 1993) and was generally maintained below 10~^. 

The value of Alr is unaffected by decay-mode-dependent variations in detector acceptance 
and efficiency, provided that the efficiency for detecting a fermion at some polar angle (with 
respect to the electron direction) is equal to the efficiency for detecting an anti-fermion at the 
same polar angle. This fact, and the high degree of polar symmetry in the SLD detector, render 
negligible. Finally, Vp was experimentally demonstrated to be consistent with zero to a 
precision of 7 x 10~^, as described above. Calculations based on polarization buildup in the 
positron DR suggested a much smaller number, Vp < 0(10"^); hence, no correction for Vp was 
applied to the data. 

Table |5| summarizes the systematic effects discussed in this section. The total uncertainties 
due to backgrounds and accelerator-induced asymmetries were much smaller than the leading 
systematic effects due to polarimetry and energy uncertainties, as can be seen by comparing the 
last three rows of Table |5[ 

Table ^ shows the run-by-run Alr results. The previously discussed Scm-dependent radiative 
correction and its uncertainty are evident in the difference between ^lr and ^lr- These five 
results show a of 7.4 for four degrees of freedom, corresponding to a probability of 11.6%. 
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Table 5: Z° event counts and corrections for all SLD run periods. Also shown are the total 
polarimetry errors (including chromaticity and IP effects) and the relative error due to the 
electroweak interference correction needed for the conversion of Alr to ^lr- Note that because 
of low statistics, a number of effects were ignored for the 1992 data. 





1992 


1993 


1994-95 


1996 


1997-98 




5,226 


27,225 


52,179 


29,016 


183,335 


Nr 


4,998 


22,167 


41,465 


22,857 


148,259 




0.0223 
±0.0099 


0.1024 
±0.0045 


0.1144 
±0.0032 


0.1187 
±0.0044 


0.1058 
±0.0017 


/bkg (%) 


1.4 
±1.4 


0.25 
±0.10 


0.11 
±0.08 


0.029 
±0.021 


0.042 
±0.032 


^bkg 




0.031 
±0.010 


0.055 
±0.021 


0.033 
±0.026 


0.023 
±0.022 


Ac (10-4) 


1.8 
±4.2 


0.38 
±0.50 


-1.9 
±0.3 


+0.03 
±0.50 


-1.3 
±0.7 


^p(10-4) 


-29 


-33 
± 1 


+24 
±10 


+29 

±43 


+28 
±69 


Ae (10-4) 




0.0044 
±0.0001 


0.0092 
±0.0002 


-0.0001 
±0.0035 


+0.0028 
±0.0014 






-1.9 


0.0 
±2.5 


2.0 
±3.0 


4.3 
±2.9 



















Vp (10-4) 


< 0.16 


< 0.16 


< 0.16 


< 0.16 


-2 
±7 


Total correction, 
AAlrMlr, (%) 




+ 0.10 
±0.08 


+ 0.2 
±0.06 


+0.02 
±0.05 


+0.16 
±0.07 


SVe/Ve (%) 


2.7 


1.7 


0.67 


0.52 


0.52 


Electroweak interference 
correction [relative (%)] 




0.3 


0.4 


0.4 


0.39 



The sin^ 9^ results are derived from the equivalence Aj^^ = Ae, which implies 

2(l-4sin^g^) 
^LR - l + (l_4sin2eeff)2- ^'^> 
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Table 6: Summary of Alr and sin^ measurements for all SLD runs (statistical and systematic 
errors are listed separately) 



Data Set 


^LR 


"^LR 


sin^ 


1992 


0.100 


0.097 


0.2378 




±0.044 ± 0.004 


±0.044 ± 0.004 


±0.0056 ± 0.0005 


1993 


0.1628 


0.1656 


0.2292 




±0.0071 ± 0.0028 


±0.0071 ± 0.0028 


±0.0009 ± 0.0004 


1994-1995 


0.1485 


0.1512 


0.23100 




±0.0042 ± 0.0010 


±0.0042 ± 0.0011 


±0.00054 ± 0.00014 


1996 


0.1559 


0.1593 


0.22996 




±0.0057 ± 0.0008 


±0.0057 ±0.0010 


±0.00073 ± 0.00013 


1997-1998 


0.1454 


0.1491 


0.23126 




±0.0024 ± 0.0008 


±0.0024 ± 0.0010 


±0.00030 ± 0.00012 


All 




0.15138 ±0.00216 


0.23097 ± 0.00027 



The average for the complete SLD data sample is 

Al^ = 0.15138 ± 0.00216, 
sin^ e'§ = 0.23097 ± 0.00027. (19) 

Small correlated systematic effects are accounted for in the calculation of this average. 
5.2.2 Leptonic Asymmetries 



The individual lepton asymmetry parameters were determined from lepton final-states (41, 42). 
The Alr measurement determined the initial-state parameter Ag, but electron polarization also 
allows direct measurement of the final-state asymmetry parameter Ai for lepton i using the 
left-right forward-backward asymmetry (Apg = ^{Ve)A£). If lepton universality is assumed, the 
results for all three lepton flavors can be combined to yield a determination of sin'^ 9^, which in 
turn can be combined with the more precise result from Alr. The event sample used for Alr was 
purely hadronic (there was a very small 0.3±0.1% admixture of t^t~ events), and hence, the left- 
right asymmetry of the lepton events constituted an independent measurement. Although the 
lepton final-state analysis described below is more sophisticated than a simple ^LR-style counting 
measurement, essentially all the information on sin^ 9^ is obtained from the left-right asymmetry 
of these events. For example in the 1996-1998 data set, the inclusion of the distributions in polar 
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Table 7: Summary of event selections, efficiency, and purity for e^e i^i for the 1997-1998 
SLD data 



Event 




Background 




Efficiency in 


# of selected 


sample 




fraction 




1 cos 6*1 < 0.9 


events 


e+e" —I- 


e+e 


0.7% T+r- 




75% 


15,675 






0.2% T+T- 




77% 


11,431 






e^e" : fi^fJ-^ : 2-7 


hadrons 


70% 


10,841 






0.9% : 2.9% : 0.9% 


: 0.6% 







angle that are essential for the extraction of the final-state asymmetries improved the resulting 
precision on sin^O^, but only to ±0.00076, compared with about ±0.00078 obtained from a 
simple left-right event count. 

The Analysis Method Table summarizes the selection efficiencies, backgrounds, and 
numbers of selected candidates for e~^e~ , n'^fx~, and t~^t~ final states. Event-by-event maximum 
likelihood fits, for each lepton flavor, were used to incorporate the contributions of all the terms in 
the polarized differential cross section and to include the effect of initial-state radiation. Photon 
exchange terms and, if the final-state leptons were electrons, t-channel contributions, were taken 
into account. All three lepton asymmetry parameters, Ag and (Ar), were obtained from 
/i^;U^ {t^t~) final states. These A^. results were combined with A^ obtained from the e~^e~ 
final state. 

Figure 12 shows the cos^ distributions for e"'"e~, and r+r^ candidates for the 1997- 

1998 data. The pre-1997 results are similar but have smaller acceptance (|cos0| < 0.8) (note 
that the SLD event totals including this older data are 22,254, 16,844, and 16,084 for the e~^e~ , 
fJ-^fJ-" and T^T~' final states respectively). 

The high- 1 cos6\ region is increasingly sensitive to the asymmetry parameters, most apparent 



in the and r+r^ distributions in Figure 12. In 1996, the upgraded vertex detector and a 

new trigger system for forward events were installed. The improved acceptance of VXD3 

allowed for efficient track finding up to |cos0| = 0.9 (|4^). The new trigger for //"'"/i" events 
covered the angular range up to |cos0| < 0.95 by requiring two back-to-back tracks that pass 
through the IP and reach the endcap WIC used for muon identification. 



40 



T/nurQniKT QTT ^^^TT T,OCQ 




200 \r 
100 [ 

Q F. I ... I ... I ... I .. . I . .. I ... I ... I ... I . 
-O.S-Oj6-Oj4-0.2 0.2 04 O.S 0.8 



cos9 

Figure 12: Distributions of the lepton cos^ for e~^e^, , and r^r^ candidates from the 

1997-1998 data set. 

Systematic Errors This measurement is statistically limited. The 1997-1998 dataset, 
which dominates the sample, typically has statistical errors between 4 and 10 times larger than 
the total systematic error for the various asymmetry parameters. The systematic errors arise 
from polarimetry, backgrounds, radiative corrections, polarization effects ("V— A"), incorrect 
charge assignment for tracks at large cos 9, and nonuniformities in the detector efficiency and 
forward-backward asymmetries. These errors are given in Table ^. 

One error deserves some discussion because of its specific relevance to polarized forward- 
backward asymmetry measurements. The dominant systematic error in the analysis resulted 
from final-state polarization effects, that introduced a selection bias. The term V— A ("vector 
minus axial vector") used in Table ^arises from the (1 — ^^)^/^ Lorentz structure of the charged- 
current Lagrangian that governs decay. For example, if both r+ and r~ decay to vri/, helicity 
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Table 8: Summary of statistical and systematic uncertainties in units of 10~^ for the 1997-1998 
data. The superscript on each asymmetry refers to the lepton sample from which it was derived 
(electrons, muons, or taus). 



Source 


At 

e 


e 


AI 

e 


Af, 


T 


Statistics 


110 


130 


130 


180 


180 


Polarization 


8 


8 


8 


8 


8 


Backgrounds 


5 




13 




14 


Radiative correction 


23 


2 


2 


3 


2 


V-A 










18 


Charge confusion 








7 


11 


Detector asymmetry 










4 


Nonuniform efficiency 


2 











conservation requires that both pions generally have lower momentum for a left-handed t~ and 
right-handed r"*" and higher momentum otherwise. This effect, which biased the reconstructed 

event mass, was significant because of the beam polarization, which induced a large and 
asymmetric polarization as a function of polar angle. The value of A^ extracted from t~^t~ 
final states was not affected, since the overall relative efficiencies for left-handed beam and right- 
handed beam events were not changed significantly (only the polar angle dependence of the 
efficiencies was changed). 

Results Results for all SLD data sets were combined, taking into account small effects due to 
correlations in systematic uncertainties (polarization and average SLC center-of-mass energy). 
Prom purely leptonic final states, one obtains Ae = 0.1544ib0.0060. This Ag result was combined 
with the left-right asymmetry measurement in the tabulation of leptonic asymmetry results: 

Ae = 0.1516 ± 0.0021, (with^^j^) 

A,, = 0.142 ± 0.015, and (20) 

Ar = 0.136 ± 0.015. 

These results are consistent with lepton universality and hence can be combined to obtain A^, 
which in the context of the standard model is simply related to the weak mixing angle. The 
result is discussed in the following section. 
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5.2.3 Combined Results for sin^ 9^ 

Assuming lepton universality, the ^lr result and the results on the leptonic left-right forward- 
backward asymmetries are combined, where small correlated systematics are accounted for, 
yielding 

= 0.15130 ± 0.00207. (21) 

This measurement is equivalent to the determination 

sin^ = 0.23098 ± 0.00026, (22) 

where the total error and corresponding systematic error of ±0.00010 are more precise than 
those obtained by any other technique. 

5.3 Quark Couplings 

It is particularly interesting to measure the couplings of individual quark flavors separately, 
and thereby probe the radiative corrections to the Zqq vertex. The b and c quarks are the 
heaviest charge 1/3 and charge 2/3 quarks, respectively, that are accessible at the energy. 
Potential new-physics signatures may preferentially appear in the heavy-quark couplings, e.g., 
those that involve Yukawa-type couplings favoring fermions with large masses. More generally, 
any unexpected difference in quark coupling of one flavor compared with other flavors could be 
a vital clue toward a solution to the puzzle of fermion family degeneracy. The relatively similar 
production rates of all quark flavors in decays, combined with our ability to tag bottom 
and charm hadron decays, offers the possibility to test the Z^ coupling to the individual quark 
flavors with high precision. In addition to the measurements with heavy quarks, comparing the 
measurements for the s quark to those of the b quark provides a unique test of the universality 
of couplings for charge 1/3 quarks. 

For most of the measurements in this section, the hadronic Z° events were selected by re- 
quiring at least seven charged tracks and a visible charged energy of at least 18GeV. The 
events were typically required to be well-contained in the high-quality tracking fiducial volume 
of I cos ^thrust I < 0.7. The flavor-tagging efficiencies referred to are generally for these selected 
fiducial events. 

5.3.1 Rb and Rc Measurements 

The formulation of Rb and Rc as ratios of hadronic cross sections ensures that the propagator 
(oblique) electroweak radiative corrections and QCD radiative corrections that are common to 
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all flavors largely cancel, isolating the heavy quark to coupling vertex radiative corrections. 
In the standard model, the large top quark mass introduces a —1.5% correction on Rf, (^3|), 
compared with the tree-level prediction. Extensions to the standard model can also produce 
potential deviations in Rf, at ~1% (^^. The effort toward < 1% precision Rf, measurements 
has therefore been one of the primary activities of the Z'^-pole experiments. The Zee vertex 
corrections in the standard model are much smaller, so that any deviations of the measured Rc 
from the standard-model prediction would signal exotic new physics processes. 

At the peak of the Rb and Rc "crisis" in early 1996, the world average for Rb was over 3a higher 
than the standard model, whereas that for Rc was over 2a lower than the standard model (p5|). 
SLD's crucial contribution was to introduce an improved analysis method that was eventually 
adopted by other experiments, resulting in a significant increase in precision. 
Rb Measurement Recent Rb measurements have generally adopted a double-tag technique 
to reduce modeling uncertainty. Events were divided into two hemispheres by the plane perpen- 
dicular to the thrust axis, and a 6-tagging algorithm was applied to each hemisphere in turn. The 
measured hemisphere tag rate Fs and event double-tag rate allow the extraction of both Rb 
and the hemisphere 6-tagging efficiency from the data by solving two simultaneous equations: 

Fg = ebRb + ^cRc + ^udsi^ — Rb — Rc), 

Fd = CbelRb + CcelRc + elasi^-Rb-Rc). (23) 
The small background tagging efficiencies for uds and charm hemispheres, e^ds and ec, as well 

^double ^double 

as the 6-tagging hemisphere correlations Cb = — and Cc = — , were estimated from the 
Monte Carlo simulation. Furthermore, a standard-model value of Rc was assumed. The Rb 
statistical error is approximately oc 1/eb, whereas the uds and charm systematic errors scale as 
Cuds/^b and ec/eb, indicating the need to maintain both high efficiency and high purity for the b 
tag. 

For the preliminary analysis of the 1996-1998 data, a cut on the neural-net c — b separation 



variable of Scb > 0.75 (see Section 4.3) was used as the b tag. Figure 13a shows the hemisphere 
tagging efficiencies and b purity for the 1997-1998 data as a function of the Scb cut. It can be seen 
that the measured 5-tagging efficiency is in good agreement with the Monte Carlo simulation. 
At the nominal cut of S'cfe=0.75, the tagging efficiencies for various flavors were euds=0.07%, 
ec=l.l%, and e^'^*'^=61.8%, which corresponded to a 6 purity of 116=98.3%. The small hemisphere 
correlation (1 — Cf, = 0.00007ib0. 00050) was due to a combination of several effects related to the 
IP, 6-tagging efficiency dependence on 6-hadron momentum, and cos 6, which were all verified 
to be small. Table M presents the preliminary Rb measurement result and systematic errors. 
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Figure 13: (a) 6-tagging efficiency and purity versus Set cut. (b) Rh statistical and systematic 
errors versus Scb cut. 

combined with the pubhshed 1993-1995 result. The choice of the Scb > 0.75 cut was the result 
of an optimization that minimized the total error, as seen in figure [T^ . 

Rc Measurement All Rb measurements now rely on the long 6-hadron lifetime to exploit 
the double-tag technique, but the shorter charm lifetime and limited vertex detector resolution 
prevent the effective use of that technique for Rc measurements in most experiments. This led to 
more complicated paths for the Rc measurements at LEP, mostly involving techniques that relied 
on fully or partially reconstructed charm hadrons (^). These methods were not only limited 
statistically, owing to low efficiency, but also limited systematically at the ±0.010 level, owing 
to the lack of clean self-calibration of tagging efficiencies from data. In 1997, SLD introduced 
an efficient double-tag analysis, which led to, by far, the most precise Rc measurement. 

The current preliminary SLD measurement of Rc is based on the 1996-1998 data only. Since 
the b tag separated out the majority of b hemispheres with high purity, the Rb analysis used a 
simple double-tag technique. The Rc measurement, on the other hand, benefited significantly 
from the additional statistics of cc hemispheres in a wider range of Scb values. Although the 
charm purity was moderate compared with the 6-tagging purity, the remaining impurities were 
mainly b hemispheres whose fractions could be precisely measured from data. The Rc measure- 
ment therefore adopted a multi-tag scheme, which was also used in some LEP Rb measurements. 
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Table 9: Hemisphere tagging efficiencies and purities for tags used in tlie Rc analysis, based on 
the c — b separation variable Scb (boxed efficiency values were measured values from data) 

Tag Name "c-pure" "c-like" "b-hke" "6-pure" 

Tag Cuts Scb < 0.3 0.3 < Scb < 0.5 0.5 < Scb < 0.75 Scb > 0.75 



eb [Vo) 
ec (%) 

^uds (%) 



2.53 



2.96 



5.10 



62.02 



17.94 



5.04 



2.29 



0.05 



0.10 



0.12 



1.12 
0.07 



b purity (%) 15.0 
c purity (%) 84.2 
uds purity (%) 0.9 



40.9 
55.3 
3.8 



70.4 
25.1 
4.5 



98.3 
1.4 
0.3 



Hemispheres were classified into five categories, corresponding to four exclusive categories of b 
and c tags defined by various Scb ranges and a fifth category with no tag. Table ^ lists the 
efficiencies and purities of the four different b and c tags for the 1997-1998 data analysis. The 
measured event rates, Nij/Ntot for two hemispheres (categories i versus j), were used in a fit 
for all tag combinations: 

AT /AT n J)^b/~ib D ^c^cfic n \ ^uds ^uds /^uds 

2 ^^ij / ^^tot - ^be^ejUij - KcejejUij - [L - Kb - ^c)^i (-^jj 

^ "V cT(iV,,/iVtot) ■ ^ ' 

Besides Rb and R^ all 6-tagging and c-tagging efficiencies (except the c-tagging efficiency for 
the 6-pure tag) were also measured from data. The small uds efficiencies and the hemisphere 
correlations (Cjj) were estimated from the Monte Carlo simulation. The hemisphere correlation 
factors were found to be very close to 1 for all tags. The essence of this multi-tag scheme is that 
by looking at the tags opposite the high-purity 6-pure and c-pure tags, the flavor composition in 
the less pure tags was measured from data. However, this efficiency cross calibration is reliable 
only when the uds fractions are kept as low as they were in this analysis. 

Table ^ summarizes the preliminary Rc result and its systematic errors. The value of the si- 
multaneously fitted Rb was also very consistent with the simple double-tag Rb result. The various 
measured b- and c-tagging efficiencies were found to be close to the Monte Carlo simulation. 



Rb and Rc Summary The SLD measurements of Rb and Rc are summarized in Table IC, 
where the results were corrected for 7-exchange by -1-0.0003 for Rb and —0.0003 for Rc. 

These results are in good agreement with the LEP measurements (48) and the standard 
model (see Figure 14). Despite the disadvantage of event statistics, the SLD R^ measurement 
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Table 10: Preliminary results and systematic summary for the Rb and Rc measurements 



Observable 



Rb 



Rc 



Data sample 
References 



1993-1998 1996-1998 



Measurement value 
Statistical error 



0.21641 
0.00092 



(|6|,|3) (im 



0.17382 
0.00308 



Monte Carlo statistics 


0.00024 


0.00095 


Event selection bias 


0.00026 


0.00027 


uds and charm physics 


0.00042 


0.00142 


6-hemisphere correlation 


0.00023 


0.00025 


g ^ cc,g -^bb 


-0.00023 


-0.00082 


Detector effects 


0.00043 


0.00079 


Rc (±0.006) 


-0.00020 




Rb (±0.0015) 




-0.00020 


Total systematic error 


0.00080 


0.00209 



is competitive in precision, and the SLD Rc measurement dominates the world average with 
a systematic error three times smaller than that of other measurements. The Rb crisis is now 
believed to have been caused by underestimated detector resolution effects that resulted in an 
excess of light-flavor tags in one of the early LEP measurements. The Rc deviation from the 
standard model was related to incorrect assumptions on charm production and decay branching 
ratios, which were later measured directly from the data. Recent measurements are less prone 
to these problems as improved analysis methods have largely reduced the sensitivity to physics 
and detector modeling uncertainties. 

5.3.2 Ab and Ac Measurements 

The direct measurements of Ab and Ac at SLD are mostly sensitive to deviations in the right- 
handed Zbb and Zee couplings, which have unambiguous standard- model predictions with small 
radiative corrections. Most proposed extensions to the standard model predict changes in the 
left-handed Zbb coupling, with little or no effect on the right-handed coupling. Any deviations 
from the standard-model prediction for Ab and Ac will therefore be sensitive to exotic new 
physics processes. 
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Figure 14: Summary of and Rc measurements from LEP and SLD. The shaded band for the 
Rf) standard model (SM) prediction reflects the uncertainty in the top quark mass. 



The heavy-quark asymmetry analyses typically involved quark flavor tags to enrich the fraction 
of the desired heavy-quark species. The primary QQ production axis was typically approximated 
by the event thrust axis. Various techniques were then used to determine the direction of the 
quark as opposed to the antiquark. We refer to this process as the quark charge assignment. 
All analyses adopted maximum likelihood fits to the polarized differential cross section (see 
Equation ^) as a function of cos0, and events were weighted appropriately based on the correct 
quark charge assignment probability, event quark- flavor fractions, and electron beam helicity. 
The information drawn from this procedure is equivalent to that obtained from the left-right 
forward-backward asymmetry at first order. The measurements were corrected for QCD 
and QED radiative effects. The dependence of radiative corrections on center-of-mass energy 
was negligible. 

In the following, we first describe the analysis procedure, emphasizing features that are special 
to SLD. The results are presented in Table |ll|. 

Af, and Ac with Leptons The use of leptons to identify heavy-hadron decays has been a 
traditional technique for the b- and c-asymmetry measurements. The lepton tags not only enrich 
the b- and c-quark events but also provide the crucial separation between quark and antiquark. 
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In this SLD measurement, the estimation of correct charge assignment probabihty for this Q/Q 
separation was based on the probabilities from the Monte Carlo for the candidate lepton to be 
from various sources, e.g., b^i, b^c^i, c^i or lepton misidentification, etc. This 
analysis benefited from the well-understood kinematics of semileptonic decays of heavy hadrons 
and the generally well- measured decay branching ratios for the various lepton sources. Besides 
the conventional lepton total momentum and transverse momentum with respect to the jet axis, 
the vertex mass b tag and other vertexing observables were also used to improve the lepton 
source classification. A special variable assisting this separation was the lepton longitudinal 
position along the vertex axis, L/D (see Section [4.2| ), which had significant resolving power 
between direct secondary 6 — > £ and cascade tertiary b ^ c ^ I leptons. The major sources 
of physics systematic uncertainty were the various semileptonic branching ratios and B mixing 
rates, which were taken from the LEP combined lepton fit results (^). 

with Jet Charge The use of momentum-weighted track charge (more commonly refered 
to as "jet charge") to sign the 6-quark direction is now a standard technique for 6- asymmetry 
measurements. The method is based on the correlation between the primary quark charge and 
the net charge of high-momentum tracks in the jet. In this analysis, bb events were selected by 
requiring at least one hemisphere with a vertex mass tag of M^orr > 2 GeV. The momentum- 
weighted track charge was calculated from 

Q= Yl qi-sign{p,-f)\{pi-f)\^, (25) 

tracks 

where qi and pi are the charge and momentum vector of track i, T is the thrust axis direction 
and K was chosen to be 0.5 to optimize the charge determination. 

The correct charge assignment probability was calibrated from data. The approximately 
Gaussian shape of the jet charge distribution allowed a simple parameterization of the correct 
charge assignment probability as Pf,(|Q|) = 1/(1 -|- e~°''l'^l). The widths of the hemisphere 
charge sum and difference distributions were used to calculate the single parameter Ob. The 
correct charge assignment probability was evaluated for each event and used in the likelihood 
fit. On average, this probability was ~69%. The small udsc background was subtracted based 
on the Monte Carlo simulation, but the b purity was measured from data using the double- 
tag technique, as in the case of the Rb measurement. There was a small hemisphere charge 
correlation, analogous to the hemisphere tag correlation in the Rb analysis, which was estimated 
from the Monte Carlo. 
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Ac with D^*^ and inclusive D*+ — > D^tt^ A traditional way of tagging charm is simply to 
reconstruct D^*^ decays. Two measurements of Ac were performed. The first analysis used fully 
or partially reconstructed D^, and mesons;^ the second used inclusively reconstructed 
D*^ D^TT~^ decays. These methods not only selected the cc events, the D^*^ i^^*^) mesons 
naturally tagged c (c) quarks with high purity. The behavior of the random combinatorial 
background (RCBG) events can be conveniently studied with the events in the mass sidebands. 
A simple b tag with mass requirement of Mcott > 2 GeV was used to veto real D^*^ background 
from bb events. 

The effective rejection of RCBG using precision vertexing allowed the inclusion of many decay 
modes in this measurement. The exclusive D*~^ D^tt^ analysis utilized four decay modes: 
i^^vr^, K^TT^TT^, K^TT^TT^TT^ , and K~£~^iJ. The direct D meson reconstructions used two 
decay modes: K^tt^tt^ and K~tt^. GRID kaon identification was used for 

the K^TT^ mode. The sample of all D^*^ candidates had the following composition: 

c-^D/b^ l)/RGBG=71%/7%/22%. 

The inclusive D*^ D^^^t analysis exploited the fact that a high-momentum leading D* 
in a cc jet would follow the jet axis very closely. Because of the low of the D*~^ — > D^-k^ 
decay, the vr^ would also travel very close to the jet axis. The tt^ candidates having momentum 
transverse to the jet axis pj^ < 0.01 (GeV/c)^ were selected with a signal-to-background ratio 
of 1:2 and an efficiency three times higher than the exclusively reconstructed sample. For the 



combined result in Table 11, the overlapping candidates between the two analyses were removed 
from the inclusive analysis for the combined Ac result. 

Af, and Ac with Vertex Charge and Kaon Tag The most precise Af, and Ac measurements 
at SLD were based on two novel quark charge assignment techniques: vertex charge and identified 
kaon charge. 

The analysis selected b and c hemispheres mainly using the neural- net c — b separation variable 
Scb (see Section 4.3). The hemisphere b tag required Scb > 0.9 and Mcorr < 7 GeV, and the c 
tag required Scb < 0.4 and a momentum sum of all secondary tracks > 5 GeV/c. 

A clean reconstruction of the secondary vertex charge tagged the heavy-quark charge. This is 
especially beneficial in the case of b hadrons because both lepton and jet-charge techniques suffer 
from dilution from neutral i?-meson mixing. This method required clean and efficient separation 
of primary and secondary tracks, fortunately achievable with the precision IP and tracking at 



SLD. The secondary track selection procedure can be found in Section Figure 15 shows the 



^Charge conjugate modes are implied throughout this paper. 
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Figure 15: Vertex charge distributions for the 6-tagged hemispheres, (o) with fully fitted tracks 
and (6) including VXD track segments. The "S^" category shown in the plot includes all neutral 
b hadrons. 

vertex charge distributions for the 6-tagged hemispheres, demonstrating good agreement between 
data and Monte Carlo and a clear +/— charge separation. Note that the inclusion of the VXD 
track segments in the charge reconstruction noticeably improved the correct charge assignment 
fraction. For the b- (c-) tagged hemispheres, the fraction of hemispheres with the secondary 
decay identified as charged was 56% (45%), and the correct b- (c-) quark charge identification 
probability was 81% (91%). 

Another quark charge assignment technique was to use the dominant 6 — c — ^ s — ^ K~ and 
c — > s — > K~ decay chains, in which the were identified by the GRID. However, in the case 
of b hemispheres, the additional contribution from the tag was found to be small. This was 
mostly due to the fact that the hemispheres with no vertex-charge tag contained an enhanced 
fraction of neutral B mesons, which had inherently poorer tag performance due to B mix- 
ing. On the other hand, a K"^ tag fraction of 30% and a correct c-quark charge identification 
probability of 86% for the c-tagged hemispheres, provided a very effective additional charge as- 
signment contribution. Therefore, the Af, analysis used the vertex-charge tag only, while the Ac 
analysis used both vertex-charge and kaon-charge tags. 

The vertex-charge and kaon-charge tags were used to make a joint quark charge assignment 
in the Ac analysis. In cases of conflict between a vertex-charge and kaon-charge assignment 
in the same hemisphere, no charge assignment was made. Thanks to the high efficiency and 
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Figure 16: Event thrust cos 9 distributions for (a) left- and (6) right-handed beams in the vertex- 
charge Ah analysis, for data (points) and Monte Carlo simulation (histograms). 

the good performance of the quark charge assignment, we were able to determine the event 
flavor composition and quark charge assignment probability simultaneously from data using a 
hemisphere double-tag technique. This procedure used the small uds efficiencies and hemisphere 
correlations from the Monte Carlo simulation and assumed the world average values of Rb and 
Rc- 

For the A), and Ac fits, events with either hemisphere having a b tag were classified as bb 
events, whereas events with either hemisphere having a c tag and no hemisphere having a b tag 
were classified as cc events. Events with two hemispheres having the same charge were discarded, 
and events with two hemispheres having the opposite charge were weighted by the joint correct 
charge assignment probability. The bb sample had a b purity of 97.5 ±0.5% and a correct 6-quark 
charge fraction of 81.5 it 0.5%. The 6-quark cos6 distributions of 6-tagged events with vertex 



charge are shown in Figure |16|. The most significant systematic uncertainty was due to charge 
assignment calibration statistics. 



Af, and Ac Summary Table 11 summarizes results for the various Aj, and Ac measurements. 
The raw measured asymmetry parameters were corrected for QCD and QED radiative effects. 
The QED corrections for Ai, and Ac were rather small: —0.2% and +0.2%, respectively. These 



were already taken into account for results in Table |ll|. The size of these corrections can be 
compared to the analogous QED corrections of 2.5% and 9% applied to unpolarized asymmetries 
measured at LEP. 
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Table 11: Summary of Af, and Ac results, detailing statistical and systematic uncertainties. All 
measurements are preliminary except for the Ac measurement with D^*\ 



Observable 



A, 



Ac 



Measurement type 
Data sample 
References 



Lepton 
1993-98 
©I) 



Jet-Q 
1993-98 



Vtx-Q 
1996-98 



(||,|2D (iiiD 



Lepton L>(*) Vtx+K 
1993-98 1993-98 1996-98 
(|0|,|5lD dH) (||,|5|) 



Measurement value 
Statistical error 



0.924 
0.030 



0.907 
0.020 



0.921 
0.018 



0.589 
0.055 



0.688 
0.035 



0.673 
0.029 



Monte Carlo statistics 
Calibration statistics 
Calibration systematic 
Detector effects 
6-tag purity 
Br(semi--£) & B mixing 
6, c decay/prod, model 
g ^ cc,g ^bb 
Beam polarization 
QCD correction 
Ac (±0.035) 
Ab (±0.035) 



0.005 



0.010 
0.009 
0.012 
0.011 
0.002 
0.005 
0.005 
0.003 



0.014 
0.016 
0.008 
0.005 
0.000 
0.000 
0.001 
0.005 
0.003 
0.001 



0.003 
0.013 
0.007 
0.003 
0.007 
0.000 
0.003 
<0.001 
0.005 
0.003 
0.001 



0.017 



0.020 
0.036 
0.022 
0.018 
0.001 
0.004 
0.005 

0.008 



0.010 
0.018 
<0.001 
0.003 
0.003 

-0.003 



0.004 
0.021 
0.009 
0.003 
0.003 
0.000 
0.004 
0.002 
0.004 
0.001 

-0.001 



Total systematic error 



0.023 



0.024 



0.018 



0.053 



0.021 



0.024 



Whereas the QED corrections were generally calculated reliably with small uncertainties, the 
QCD corrections had significant theoretical uncertainties and were sensitive to analysis details. 
The total integrated QCD correction was at the ~ 3.5% level for an ideal measurement of 
either A^^ or Aq. However, the analysis event selection and final-state particle reconstructions 
typically tended to suppress the events with hard gluon radiation and therefore only needed 
reduced QCD corrections. This can be understood because b tags were typically more efficient for 
high-momentum b hadrons, and high-x D* reconstruction and high- momentum lepton selections 
also favored high-momentum c and b quarks in two-jet events. The QCD corrections were 
scaled down by analysis-dependent factors ranging between 25% and 75%. The SLD A^ and 
Ac analyses used a first-order QCD correction calculated by Stav & Olsen (|56|), which included 
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heavy-quark mass effects with exphcit cos 9 dependence. The second-order QCD corrections and 
their uncertainties were based on the prescription found in Reference (|57|). 

The individual and Ac measurements were combined, taking into account systematic cor- 
relations. The various A^, measurements had significant event sample overlaps. A statistical 
correlation matrix was built, taking into account the different weight of each event used in the 
analysis. The effective correlations obtained from this procedure were as follows: 

Lepton versus Jet-Q: 22%, 
Lepton versus Vtx-Q: 15%, 
Jet-Q versus Vtx-Q: 32%. 

The combined result also included an older measurement of Ai, using the -fC tag alone from 



the 1994-1995 data (55). Given the much smaller samples for the various Ac analyses, their 
statistical correlations were deemed to be small. The combined preliminary SLD At, and Ac 
results are 

Ah = 0.916 ±0.021, 

Ac = 0.670 ±0.027. (26) 

Figure ^ lists the individual Ai, and Ac measurements, as well as the SLD average, along with 
the indirect measurements derived from the LEP Aps measurements, assuming a measured A^ 
from the SLD and LEP combined ^lepton result of Ac = 0.1501 ± 0.0016. 

5.3.3 As Measurement 

An important test of the standard model is to verify the universality of coupling strengths of 
the boson to fermions with the same charge and weak isospin. As described in the previous 
section, a number of precise measurements of heavy-quark asymmetries have been performed, but 
few measurements of light-quark asymmetries have been carried out, owing to the lack of clear 
experimental signatures to separate uu, dd, ss decays from one another. Nevertheless, it is 

possible, albeit with low efficiency, to select a fairly pure sample of Z^ — > ss decays by requiring 
each event hemisphere to contain high- momentum kaons, as outlined below. 

In the As analysis (^), Z^ — s- hb and Z^ — s- cc decays were suppressed by requiring the 
events to contain no more than one well-measured track with normalized impact parameter in 
the transverse plane d/ad > 2.5. Charged kaons were identified with the GRID and required to 
have p > 9 GeV/c to provide a sample 91.5% pure in K^. Neutral kaons were identified via their 
decay into tt'^tt^ final states with a purity of 90.7%. An event was tagged as Z^ — > ss if one 
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Measurements 



Measurements 



SLD JetC 
SLD Lepton 
SLD K* tag 
SLD VtxQ 
SLD Average 

ALEPH Lept 
DELPHI Lept 
L3 Lept 
OPAL Lept 
ALEPH JetC 
DELPHI JetC 
L3 JetC 
OPAL JetC 
DELPHI NN 
LEP Average 







!+-•- 


H 

1 






H • 1 




<^ SM 



0.907 ±0.020 ±0.024 


SLD soft Jt' 


0.924 ± 0.030 ± 0.023 


SLD d',D* 


A OCC M A AOO -1- A H AO 


SLD Lepton 


0.921 ±0.018 ±0.018 






SLD K S vtx-0 


0.916 ±0.021 






SLD Average 


0.886 ±0.035 ±0.020 




0.918 ±0.052 ±0.022 


ALEPH Lepton 


0.873 ±0.058 ±0.029 


DELPHI Lepton 


0.851 ±0.038 ±0.021 


L3 Lepton 


0.911 ±0.024 ±0.014 


OPAL Lepton 


0.892 ±0.042 ±0.01 6 




ALEPH D* 


0.843 ±0.090 ±0.050 




0.894 ±0.049 ±0.036 


DELPHI D' 


0.883 ±0.032 ±0.021 


OPAL D' 


0.880 ±0.018 


LEP Average 



-SM 



0.685 ± 0.052 ± 0.038 
0.690 ±0.042 ±0.021 
0.589 ± 0.055 ± 0.053 
0.673 ± 0.029 ± 0.024 

0.670 ± 0.027 

0.580 ± 0.047 ± 0.040 
0.645 ±0.080 ±0.061 
0.774 ±0.314 ±0.160 
0.575 ± 0.054 ± 0.039 
0.61 7 ±0.080 ±0.024 
0.635 ±0.083 ±0.025 
0.628 ±0.104 ±0.050 

0.608 ±0.032 
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Figure 17: Comparison between the A^, and Ac measurements by SLD and the A), and Ac values 
derived from the indirect Ap^ and Ap^ measurements at LEP. 

hemisphere contained a candidate and the other contained either an oppositely charged 
or a Kg candidate. The respective ss purities for K^K^ and K^K^ events were estimated to 
be 73% and 60% from the Monte Carlo simulation, with an overall ss event selection efficiency 
of 2.6%. Identification of the s-quark hemisphere relied on the charge and was correctly 
assigned with a probability of 97.5% for K~^K~ events and 85.0% for K^K^ events. 

A maximum likelihood method was used to extract Ag from the s-quark cos 9 distributions for 
left- and right-handed electrons in both samples (see Figure ^). From the figure, it is apparent 
that bb and cc backgrounds have asymmetries with the same sign as the signal asymmetry, 
but that uu and dd events contribute with the opposite sign. Therefore, uncertainties in the 
modeling of light quarks have the greatest impact on the value of Ag. Using the 1993-1998 data, 
the analysis determined Ag = 0.895 it 0.066 (stat.) it 0.062 (syst.), where the systematic error is 
dominated by uncertainties in the uu and dd contributions to the cos 9 distributions. To reduce 
modeling uncertainties, these contributions were constrained using the data, as was the s-quark 
correct charge assignment probability. 
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Figure 18: Distributions of the s-quark cos 9 for data (points) in (a) K events with Ve < 0, 
(b) K+ K- events with Ve > 0, (c) K° events with Ve < 0, (d) events with Ve > 0. 

The histograms represent the result of the fit. 

6 b^-b2 mixing 



6. 1 Introduction 

The special features of the SLC and SLD environment were particularly well-suited for studies 
of B decays. Measurements of the average 6-hadron lifetime (|59| ) and of the and B^ life- 
times (^) provided the initial impetus for the development of the inclusive topological vertexing 
described above. Recent updates of the B^ and B^ lifetime measurements (|6l| ) are among the 
most precise available to date and have helped establish a small but significant difference between 
the lifetimes of the B~^ and B^ mesons, in agreement with predictions based on the heavy-quark 



expansion technique (|62|,|63D. 

Investigations of the properties of B meson decays at SLD have included the search for charm- 
less hadronic decays into two-body and quasi-two-body final states (]6^), the search for flavor- 
changing neutral-current 5 — > s gluon decays (^), and the measurement of the charm hadron 
yield ®. 



A topic of particular interest is the study of the time dependence of B^-B^^ and B^-B^ mixing. 
Experimental studies of mixing exploit all the special features of the SLC/SLD environment, i.e., 
beam polarization, small beam spot, high-resolution vertexing, and good particle identification. 
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Below, we focus on B'^-B^ mixing studies in some detail. A recent study of B'^-B^ mixing with 



a kaon tag is described elsewhere (|67|). 
6.2 Theory 

The neutral B meson system consists of B^ and B^ flavor eigenstates, which are superpositions 
of heavy and light mass eigenstates {q = d and s for B'^ and i?^ mesons, respectively). The 
mass eigenstates evolve differently as a function of time, resulting in time-dependent B^j-B^ 
oscillations with a frequency equal to the mass difference Aniq between the heavy and light 
eigenstates. As a consequence, an initially pure \Bg) state may be found to decay as \Bg) or 
\B^) at a later time t with a probability equal to 

r 
— ( 

2 



P{BI^bI) = -e-^*(l + cosAm,t) or (27) 



P(SO^i?0) = ^e-r*(l-cosAm,f). (28) 



(Here we neglected the lifetime difference between mass eigenstates. 
The oscillation frequency can be computed to be (|68|) 



Am, = ^mBm'iF{mi/M^)flBB,m \V:,Vt,\^ , (29) 

where Gp is the Fermi constant, ms, is the B^ hadron mass, nit is the top quark mass, My/ 
is the W boson mass, F is a function defined in Reference (|69|), and r/qcD is a perturbative 
QCD parameter. The parameter B^^ and the decay constant f^q parameterize hadronic matrix 
elements. For further details, see the review by Gay (f70|). Much of the interest in B mixing stems 
from the fact that a measurement of the B^ {B^) oscillation frequency allows the magnitude of 
the poorly known CKM matrix element Vtd (Vts) to be determined (see Equation However, 
the extraction of \Vtd\ from the fairly precisely measured value Ama = 0.472 it 0.017 ps~^ (^) is 



affected by a theoretical uncertainty of about 20% in the product /b, y -Bs, (72). Uncertainties 
are reduced for the ratio 



Arus msJlBB, 



V, 



ts 



Vtd 



2 



-(1.16 ±0.05)^ 



Vt 



ts 



Vtd 



2 



(30) 



which indicates that the ratio |V(s/Vtrf| can be determined with an uncertainty as small as 
5%(|7D. 

The CKM matrix describes weak quark mixing in the standard model and contains an irre- 
ducible phase that provides the mechanism for CP violation in weak decays. This matrix can be 
written in terms of four fundamental parameters: A = sin^cabibboi A, p, and rj [as defined in the 



SLD PHYSICS 



57 



Wolfenstein parameterization ([73|)]. The value of ry is directly related to the amount of CP viola- 



tion. The parameters A and A are well-measured, A = 0.2237±0.0033 and A = 0.819±0.040 (|7|), 
but p and ij are among the least-known parameters of the standard model. In terms of these 
parameters, we have Am^ oc \Vtd\'^ — A^A^[(1 — p)^ -|- r/^] and Arris oc |Vfsp ~ j4^A^, where 
p = p{l — ^) and f] = r]{l — As a result, the combination of and Bg mixing provides one 
of the strongest constraints on p and r], and therefore on CP violation, in the standard model. 

6. 3 Experimental Ingredients 

Experimental studies of B^-B^ mixing require two main ingredients: (a) the measurement of the 
Bg decay proper time, and (b) the determination of the B^ or B^ flavor at both production and 
decay to classify the decay as either "mixed" (if the tags disagree) or "unmixed" (otherwise). 
The significance for a B^ oscillation signal can be approximated by ( [75|) 



where N is the total number of decays selected, fs is the fraction of B^ mesons in the selected 
sample, w is the probability of incorrectly tagging a decay as mixed or unmixed (i.e., the mistag 
rate), and at is the proper time resolution. The proper time resolution depends on both the decay 
length resolution and the momentum resolution Up according to = [ai/'^fic)'^ + (tap/p)'^. 
Based on the Wolfenstein parameterization, we see that Arris/ Arrid — 1/A^, which is of the 
order of 20. Therefore, B'^ oscillations are expected to be much more rapid than i?^ oscillations. 
The ability to resolve such rapid oscillations thus requires excellent decay length resolution and 
benefits from having a low mistag rate and a high B^ purity. The importance of decay length 
resolution is illustrated in Figure |l9|, which shows the fraction of 6-hadron decays tagged "mixed" 
as a function of proper time for Arris = 20 ps~^, B^ purity of 18%, mistag rate w = 0.25, relative 
momentum resolution ap/p = 10%, and either ai = 200 /um or ai = 60 pm, the former being 
typical for LEP experiments and the latter typical for SLD. It is clear from Figure ^ that the 
decay length resolution is the most critical ingredient for sensitivity to the large Bg oscillation 
frequencies that are expected. 

6.4 Analysis Methods 

SLD performed three separate analyses using the 1996-1998 data: "Z?s-|-tracks," "lepton-!-!)," 
and "charge dipole." These analyses differ in the way the Bg decay candidates were reconstructed 
and the flavor at decay was tagged (see below). 




(31) 
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Figure 19: Fraction of 6-hadron decays tagged "mixed" as a function of proper time for two 
different decay length resolutions, = 200 /im and 60 /xm. Other relevant parameters are 
Arus = 20ps"\ fs = 18%, w = 0.25, and ap/p = 10%. 

Common to all analyses was the determination of the production flavor (also known as the 
initial-state tag). This was achieved by combining several tagging techniques. The most powerful 
technique relied on the large polarized forward-backward asymmetry in — > bb decays and 
was unique to SLD. In this case, a left- (right-) handed incident electron tagged the forward 
hemisphere quark as a 6 (6) quark. The probability of correctly tagging a b quark at production 
is expressed as 

/ ^^ 1 , — Ve cos 9 
P,(cos(l) = - + A^-^^j—^^. (32) 

The average mistag rate was 28% for an average electron beam polarization of 0.73. Other 
tags were also used. They relied on charge information from the hemisphere opposite that of 
the Bg decay candidate (i.e., the hemisphere expected to contain the other b hadron in the 
event): (a) momentum- weighted jet charge, (6) secondary vertex charge, (c) charge of kaon 
from the dominant decay transition 6 — > c ^ s, (d) charge of lepton from the direct transition 
b ^ i~ , and (e) charge dipole. The first four of these tags were used to measure Af, and/or 
Ac, as described in Section |5.3| , and the charge dipole tag is described below. These tags were 
estimated to have efhciencies (mistag rates) of 100% (66%), 43% (75%), 16% (74%), 9% (74%), 
and 17% (70%), respectively. The five tags were combined to yield an efficiency of 100% and a 
mistag rate of 28%. Overall, the average mistag rate was estimated to be 22%-25%, depending 
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Initial State b-quark Probability 

Figure 20: Distribution of the computed initial state 6-quark probability for data (points) and 
Monte Carlo (histograms), for the events selected in the charge dipole analysis. 



on the cos 9 acceptance of each particular analysis. Figure 20 shows the combined 6-quark 
probability distributions for data and Monte Carlo in the charge dipole analysis and indicates a 
clear separation between b and b quarks. 

Also common to all analyses was the method to estimate the 5-hadron momentum. It utilized 
the following information: the momentum of all secondary tracks, the energy deposited in the 
electromagnetic section of the calorimeter, and kinematical constraints stemming from the known 
6-hadron mass and the e~^e~ collision center-of-mass energy. 



6.4.1 D^+Tracks Analysis 



The Z)s+tracks analysis (76) aimed to partially reconstruct the decay — > D~X with full 
reconstruction of the cascade decay D~ — > (j)Tr~ or K*^K~ . As a result, the analysis achieved a 
high purity and excellent decay length resolution but at the price of fairly low efficiency. Dg 
candidates were selected with a neural-network algorithm combining kinematical quantities and 
particle identification information to yield a sample of 280 (p-K" and 81 — > K*^K~ 

candidates. These Dg candidates were vertexed with one or more charged tracks to form the B^ 
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decay. A (D+) tagged the decay flavor as 5° (^O). Table [l| summarizes the performance 
of this analysis. Increased sensitivity to decays was obtained by separately parameterizing 
decays with a reconstructed track charge sum equal to or different from zero. For example, 
B'^ purities above 50% were obtained for the neutral sample. Furthermore, B^ — > D^l'^uX 
candidates were identified and assigned both higher purity and lower decay flavor mistag rate 
than the average. 

6.4.2 Lepton+Z? Analysis 



The lepton+Z) analysis (77) attempted to reconstruct B'^ Djt^v decays, where the D~ decay 
was not reconstructed exclusively but rather with the inclusive topological vertexing technique 
described in earlier sections. The B^ vertex was reconstructed by intersecting the lepton trajec- 
tory with the D trajectory, which is estimated from the inclusively reconstructed D vertex and 
its net momentum vector. Direct leptons from b ^ i transitions were selected with a neural- 
network algorithm that relied on the following variables: the transverse momentum of the lepton 
with respect to the B vertex direction (vector stretching from the IP to the B vertex), the B 
decay length, the transverse momentum of the lepton with respect to the D vertex direction 
(vector stretching from the B vertex to the D vertex) , the mass of the charged tracks associated 
with the B decay, and the distance of closest approach of the lepton to the B vertex. This 
algorithm effectively suppressed leptons from b ^ c ^ i transitions, which produce leptons with 
sign opposite that of direct leptons and thus dilute the decay flavor tag. To enhance the B^ 
purity, the sum of lepton and D vertex track charges was required to be zero. The purity was 
further enhanced from 16% to 39% in the subsample containing an opposite-sign lepton-kaon 
pair. Table || gives performance parameters. 

6.4.3 Charge Dipole Analysis 



The charge dipole analysis (|77| ) was based on a novel method to inclusively reconstruct decays 
of the type B^ DjX. For each candidate hemisphere, events were selected in which both 
a secondary (B) and a tertiary (D) vertex could be resolved. A "charge dipole" 5Q was then 
defined as the distance between secondary and tertiary vertices signed by the charge difference 
between them such that 6Q < {6Q > 0) tagged B^ (B^) decays (see Figure 21). This analysis. 



first developed by SLD, relied heavily on the ability to recognize secondary tracks (from the 
-B-decay point) versus tertiary tracks (from the -D-decay point). It also benefited from the fact 
that B'^ decays yield mostly charged secondary and tertiary vertices, whereas B^ mesons tend 
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X Q = +1 orO 



Figure 21: Sketch illustrating the decay cascade structure for X decays. 

to produce neutral vertices since they decay mostly into neutral D mesons. 

As in the lepton+D analysis, the total charge of all secondary and tertiary tracks was required 



to be zero, in order to enhance the purity from 10% to 15%. Figure 22 shows the 5Q 
distribution, illustrating the clear separation between B and B decays. The average decay 
flavor mistag rate was estimated to be 24%, driven mostly by decays of the type B^ DDX 
(representing approximately 20% of the selected sample) for which the method assigned the 
decay flavor randomly. In contrast, the mistag rate was estimated to be 12% for decays of the 
type B^ DjX, where X does not include any charm hadron. 

In all the analyses, the sensitivity to B^ oscillations was enhanced by evaluating the production 




-0.75 -0.5 -0.25 0.25 
Charge Dipole (cm) 

Figure 22: Distribution of the charge dipole for data (points) and Monte Carlo (solid histogram). 
Also shown are the contributions from h hadrons containing a b quark (dashed histogram) or a 
b quark (dotted histogram). 
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Table 12: Summary of relevant parameters for the mixing analyses 





1)5.+ Tracks 


Lepton+D 


Charge Dipole 


No. decays 


361 


2087 


8556 


Bs purity 


38% 


16% 


15% 


Production flavor mistag rate 


25% 


25% 


22% 


Decay flavor mistag rate 


10% 


4% 


24% 


(Jl (60% core Gaussian) 


48 /im 


55 ^m 


76 ^m 


(Tl (40% tail Gaussian) 


152 ^m 


217 urn 


311 /im 


(7p/p (60% core Gaussian) 


0.08 


0.06 


0.07 


(7p/p (40% tail Gaussian) 


0.19 


0.18 


0.21 


Sensitivity 


1.4 ps^i 


6.3 ps~^ 


6.9 ps~i 


aA at Ams= 15 ps""*^ 


1.60 


1.14 


1.11 



and decay flavor mistag rates, and the decay length and momentum resolutions, an event-by- 
event basis. 

6. 5 Results 

The study of the time dependence of B^-B^ mixing utilized the amplitude method (|75|). Instead 
of fitting for Anig directly, the analysis was done at fixed values of Arris and a fit to the amplitude 
A of the oscillation was performed, i.e., in Equations |2^ and ^8], for the unmixed and mixed 
probabilities, the term [1 it cos(Amst)] was replaced with [1 ziz Acos{Amst)]. This method is 
similar to Fourier transform analysis and has the advantage of facilitating the combination 
of results from different analyses and different experiments. One expects the amplitude to 
be consistent with A = for values of the frequency sufficiently far from the true oscillation 
frequency, whereas the amplitude should peak at j4 = 1 at that frequency. 

The L's+tracks, lepton+D, and charge dipole analyses were combined, taking into account 
correlated systematic errors. Events shared by two or more analyses were only kept in the 
analysis with the highest sensitivity, thereby removing any statistical correlation between the 
analyses. The dominant contributions to the systematic error were found to be due to the 
uncertainties in the B^ production rate, and the decay length and momentum resolutions, as 



well as the decay flavor tag in the case of the charge dipole analysis. Figure 23 shows the measured 



amplitude as a function of Arug for the combination. The measured values are consistent with 
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Figure 23: Measured amplitude as a function of Am^ for the L)c,+tracks, lepton+D, and charge 
dipole analyses combined. 

A = for the whole range of Am^ up to 25 ps~^, and no significant evidence is found for a 
preferred value of the oscillation frequency. The following ranges of oscillation frequencies 
are excluded at 95% CL: Arris < 7.6 ps~^ and 11.8 < Am^ < 14.8 ps~^, i.e., the condition 
A + 1.645 cTyi < 1 is satisfied for those values. The combined sensitivity to set a 95% CL lower 
limit is found to be at a Am^ value of 13.0 ps^^, which corresponds to the frequency that 
satisfies 1.645 cjyi = 1. All results are preliminary. 

The impact of the SLD analyses on the search for oscillations is apparent in Figure 24, 
which shows the uncertainty in the measured oscillation amplitude as a function of Arrig for 
CDF, LEP, SLD, and the world average. The growth of a a with increasing Am^ is due to the 
finite proper time resolution and clearly illustrates that the excellent proper time resolution for 
SLD yields a smaller increase in uncertainty than in the case of LEP, for example. The SLD 
sensitivity of 13.0 ps~^ can be compared to the combined LEP sensitivity of 14.5 ps~^, obtained 
with a data sample 30 times larger than that used by SLD. 
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Figure 24: Amplitude uncertainty as a function of Am^ for CDF, LEP, SLD, and the world 
average. The dashed horizontal line corresponds to the condition 1.645 = 1 used to define 
the 95% CL sensitivity. 

7 INTERPRETATION OF RESULTS 



7.1 B^^-B^^ Mixing 

As mentioned above, the study of oscillations is motivated by the fact that the oscillation 
frequencies Am^ and Am^ provide very powerful constraints on the weak quark mixing sector 
of the standard model in general, and on CP violation in particular. These constraints can be 
represented as circular bands centered around the point (1,0) in the p-i) plane (see Figure p5|) . 
Also shown are the constraints from measurements of the CP-violating parameter eK in the 
K^-K^ system and the ratio of CKM elements |Ki6/K:6! measured in charmless B decays. The 
combination of these constraints provides a measurement of the parameters p and rj. These 
parameters correspond to the apex of the Unitarity Triangle (shown in Figure 25). The triangle 
represents one of the unitarity conditions that the CKM matrix must satisfy. A nonvanishing 
value of r] implies the existence of CP violation, whereas the angles of the triangle are related 
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Figure 25: Constraints on the parameters p and rj from the measurements of ex, \Vub/Vcb\, Am^, 



and Anis- Figure extracted from Ciuchini et al (74) 



to the single CP-violating phase of the CKM matrix. 

A global fit within the context of the standard model yields p = 0.224 it 0.038 and fj = 
0.317 lb 0.040 ([74|). The fit uses the entire oscillation amplitude spectrum rather than just 
the lower limit on Am^ to incorporate information about the exclusion significance as a function 
of Arris. It should be noted that the treatment of errors in such a fit remains controversial 
because theoretical uncertainties dominate. Of special interest are the predictions for the CP- 
violating angles of the Unitarity Triangle a, (3, and 7, which are currently the subject of intense 
experimental activity. The fit yields sin 2/3 = 0.698 it 0.066, sin 2a = —0.42 it 0.23, and 7 = 
(54.8 lb 6.2)° ([74|). CP violation in B decays has now been unequivocally observed by the BaBar 
and Belle collaborations, with a world average value for sin 2/3 of 0.79 it 0.10 (|7^), in good 
agreement with the above standard-model prediction. In the next few years, it is likely that we 
will learn the most from the combination of sin 2/3 and Arris measurements at the B factories 
and the Tevatron, since they provide orthogonal constraints on the Unitarity Triangle and are 
fairly clean theoretically and experimentally. 



7.2 Electroweak - Quark Couplings 

SLD measurements of R^, Rc, Af,, Ac, and Ag are all in good agreement with the standard 
model. The SLD R^ measurement precision of ~0.5% is in the interesting regime for testing the 
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Table 13: Left- and right-handed Zbb and Zee couphngs determined from a fit to SLD and LEP 
measurements; the standard model (SM) prediction is also provided 



Coupling constant Fit result SM prediction 

g'l -0.4183 ±0.0015 -0.4211 

5^ -0.0962 ± 0.0064 -0.0774 

gl 0.3443 ± 0.0037 0.3467 

gji 0.1600 ±0.0048 0.1548 

~1% level effect due to physics beyond the standard model through radiative corrections. The 
As measurement, showing good agreement with the Af, measurement, confirms quark coupling 
universality at the ±10% level. These measurements are also mostly in agreement with the 
measurements from LEP. In the case of A), and Ac, the indirect measurements from LEP derived 
from -/Ip-g SjIicI -^pg 

are both lower than the SLD average at the 1.3cj and l.Scj levels, respectively. 
However, the derived LEP Aj, using the LEP and SLD combined A^ for the initial state A^, is by 
itself 3.1cr away from the standard model. This is essentially the same effect as the case of Ap-g 
and Apg measurements giving higher sin^ 9^ values than other measurements when assuming 
universality in the fermions (see next subsection). 

An interesting perspective on the interconnection between the measurements of Ai, Ab, and 
Apg is the Zbb vertex coupling analysis scheme from Takeuchi et al (^). Deviations {5g\,5g^^ 
from the standard model Zbb coupling can be transformed into a pair of variables {5,b,Cb) such 
that Rb deviations depend only on and Ab deviations depend only on Qb- Particularly inter- 
esting is the relation between the parity violation type variable Cfe and 5s\v?9^, as shown in 
Figure The consistency between the various measurements and the standard model is only 
at the 1.0% CL. The projection of the different bands on the horizontal line corresponding to 
Qb = 0, i.e., assuming the standard-model Zbb coupling, illustrates the discrepancy in svc? 9^ 
between the values derived from Ai and from Apg. The projection along the vertical line rep- 
resenting the measured siv? 9'^ from Ai illustrates the difference in Ab, comparing the direct 
measurement from SLD and the indirect measurement from ^pg. 

A general fit for the left- and right-handed Zbb and Zee couplings was performed (^) in the 
context of the standard model (with sin^ 9^ extracted from the world average Ai) using the 



Rb, Rc, Ab, Ac, Ap^, and Apg measurements. Table 13 summarizes the results of the fit and 



the corresponding standard model predictions. These results are also illustrated in Figure 27 
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Figure 26: The Zbb coupling analysis result following Takeuchi et al. The standard-model 
point at (0,0) is defined by mf=174GeV, mj/=300GeV, as=0.119, and a=l/128.905. The thin 
horizontal band around (0,0) corresponds to the standard model rrit and niH variations indicated 
in the plot. The error ellipses represent 68% and 95% CL contours for the fit. 



There is generally good agreement between the measurements and the standard model for the 
Zee coupling at the current precision level. The apparent departure of the Zbb coupling from 
the standard model, driven primarily by the LEP Ap^ measurements, is mainly affecting the 
right-handed coupling value. This is particularly difficult to accommodate because no known 
model can produce a deviation from the standard model at this level. The direct measurement 
of Af, from SLD provides an important constraint on a possible Zbb coupling anomaly, which has 
relevance to the validity of setting constraints on the Higgs mass from measurements of ^fb; as 
is discussed next. 
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Figure 27: Left- and right-handed Zbb and Zee couphngs combining measurements from SLD 
and LEP. The standard model (SM) predictions are also shown with small arrows indicating the 
effect of the uncertainty in the top quark mass. 

7.3 Electroweak - Lepton Couplings 

The SLD measurement of sin^ d^, precise and statistics-dominated, represents a benchmark for 
determinations of the weak mixing angle. What do we learn from this result? We must first 
examine its consistency with the standard model. We have discussed the role of radiative effects 
and noted that the measured top mass agrees with theoretical constraints at the level of about 
7 lb 10 GeV (a relative 4 it 6%). Now only one standard-model parameter remains undetermined 
(although direct experimental constraints exist) — the Higgs boson mass rriH- Because of radia- 
tive effects on the gauge boson propagators, corrections due to the Higgs boson are expected to 
contribute logarithmically in uih, and hence sensitivity to this parameter is marginal. 

7.3.1 Higgs Mass Constraints 

Nevertheless, the SLD measurement error is small enough that a meaningful constraint is avail- 
able from this measurement alone. The sensitivity of an observable O to mn can be quantified by 
the normalized partial derivative giog-^^Ql'^^ /gcY) ' '^^^re 5n is the total uncertainty in Q. The 
value of this metric for the SLD sin'^ 6^ measurement, at 4.4, is larger than any other presently 
available electroweak result [closest are the forward-backward b quark asymmetry results from 
LEP (3.9) and the W boson mass result from LEP-H (3.7)]. Figure p^ a illustrates the depen- 
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Figure 28: (a) Leading-order radiative effects on sin^ 9^ as a function of the Higgs boson mass. 
(b) Change in curve for a fit to the Higgs mass using the SLD sin^ 0^ result. 



dence of sin 9^ on the Higgs mass. From this figure it is apparent that the SLD sin 9^ result 
of 0.23097 prefers a low Higgs mass, and that mass constraints set in this region benefit from the 
steeper slope of the curve. The largest uncertainty in these standard-model calculations is due 
to the present top quark mass error of 5.1 GeV, equivalent to Jsin^ 9^ ~ ±0.00016. In addition, 
an ambiguity in the determination of the standard-model Higgs boson mass arises from uncer- 
tainties in the evolution of the fine structure constant from low to its value at the scale, 
a(M^). This uncertainty stems from the need to use low-energy e~^e~ annihilation cross-section 
data in the calculation of hadronic loop effects. Recent data from the Beijing Electron Syn- 
chrotron (BES) e~^e~ experiment (p|) have led to improved errors, which presently correspond 
to an uncertainty in sin^ 9^^ of about ±0.0001 (|82|) . However, a number of evaluations of a(M^) 
are in common use, leading to variations in sin^ 9^ predictions over a range of ~ 0.00015, and 
an associated range in quoted Higgs boson mass limits. The choice of a(M^) used throughout 



this review is due to Burkhardt and Pietrzyk (82). 

We perform a fit to the standard model by including at least one electroweak observable, 
along with the fit parameters: mn, Mz measured at LEP, nit measured at Fermilab, a{M2;), 
and the strong coupling constant Og (the Fermi constant Gi;' is held fixed). With the SLD 
sin^ 9^ result as the sole electroweak observable included in the fit. Figure p8^ gives the results 
in the form of a Higgs mass x^- The resulting one-sided confidence upper limits from the fit are 
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ruH < 133 GeV (95% CL) and < 205 GeV (99% CL). Also shown is the current direct search 
lower limit from LEP-II of 114.1 GeV. These direct and indirect Higgs bounds are in modest 
agreement: the value of the one-sided confidence limit at 114.1 GeV corresponds to a probability 
of 7%. 

As we have seen, the SLD data are consistent with lepton universality. Measurements at 
LEP of Ap^, and of polarization (which separately provides and A-r), confirm lepton 
universality as well. All available A£ data, expressed in terms of sin^O^, are consistent (see 
Figure 29). Also shown are the LEP sin"^ results deriving from quark asymmetries, which, 
taken together, provide an average (0.23230 it 0.00029) that is 3.3cj different from the lepton 
asymmetries average (0.23113 it 0.00021). Unresolved issues in the Zbb coupling data were 
discussed in the previous section. The for the overall sin^ 9^ combination has a probability 
of 2.5%. The quark asymmetry average is dominated by the Ap^ measurement, the most precise 
sin^ 9^ determination at LEP. 

Fits are made to worldwide electroweak results (^), including all SLD and LEP measurements, 
LEP-II and Tevatron Mw results, deep inelastic neutrino scattering data, and measurements of 
atomic parity violation in cesium and thallium (^). Among these, the most sensitive to itih 
are, in order of decreasing sensitivity, the asymmetries, M^y, and the total width Tz- With 
the exception of the data mentioned above, all of these measurements agree satisfactorily 



and favor a light Higgs boson. Figure 30 shows the individual contributions to the overall fit 
broken down into SLD/LEP leptonic and LEP hadronic contributions to sin^ 9^^, the world M\y 
average, and all other electroweak results (dominated by Tz)- Including all data results in limits 
of rriH < 195 GeV (95% CL) and < 260 GeV (99% CL), in agreement with the direct lower Higgs 
mass limit at 114.1 GeV. The psr degree of freedom for this fit is 23.5/16, corresponding to 
a probability of 10%. Excluding the A^^ result from the fit increases the probability to 47% 
and substantially tightens the Higgs mass limits to mn < 133 GeV (95% CL) and < 185 GeV 
(99% CL). The confidence level of this restricted fit at the value of the direct lower Higgs mass 
limit is 9%. Within the context of the standard model, the electroweak data are suggesting that 
the Higgs boson is relatively light and is in the region expected by supersymmetric extensions 
of the standard model (^). Although these constraints are not model-independent, they do 
significantly restrict possible new physics, even if the Higgs mass limits are evaded (^5|). 
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Figure 29: SLD and LEP sin^ 6^ results. The first three results are based on lepton asymme- 
tries including the SLD result, the last three results on 6-quark, c-quark, and hadronic-event 
asymmetries are based on measurements at LEP. 



7.3.2 Extensions of the Standard Model: The S, T, and U Parameters 

It is desirable to consider the electroweak data in a more general context than the standard 
model. The framework of our analysis will follow Lynn et al. (|86[), who made the following 
assumptions: 

• The gauge group of electroweak physics is SU{2)l x C/(1) {Z,Z' mixing, for example, is 
not incorporated). 

• The symmetry-breaking sector has a global SU{2) "custodial" symmetry that insures p = 
M'^/Ml cos2 ew = l + 0{a). 

• Vacuum polarization effects dominate compared to vertex and box corrections (the oblique 
hypothesis). 
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Figure 30: Change in as a function of the Higgs mass for various measurements separately. 
"Other" refers here to ah the other electroweak measurements included in our global fit. 



• Any new particles are heavy, so that Taylor expansions are appropriate descriptions of the 
vacuum polarization corrections. 

With these assumptions, the total number of independent gauge boson correction terms is six. 
Using the three precision constants {Gf, Mz, and a), this leaves three undetermined parameters. 



which are denoted S, T, and U in the scheme of Peskin & Takeuchi (87). The S parameter is a 
measure of the weak isospin-conserving new heavy sector. The T parameter is equal to a^^Ap 
and is a measure of weak isospin breaking in the heavy sector (for example, the top quark to 
bottom quark mass splitting). The U parameter introduces small effects, is only relevant for the 
observable Myy, and is ignored here. A heavy fermion doublet changes the T parameter, which 
is quadratic in the mass splitting. The standard-model Higgs boson contributes offsets that are 
logarithmic in the Higgs mass for both S and T, though these effects are opposite in sign. 

We performed a fit to the electroweak data, where offsets from S and T = were relative to a 
reference value of the top quark and Higgs masses, which we chose to be 175 GeV and 100 GeV, 



respectively. Figure 31 gives the result. 
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Figure 31: Global electroweak fit in the S and T plane. The elhpses correspond to the 68% and 
90% CL contours of the global fit. Four of the contributors to the fit are shown as ±lo" bands. 



The fit ellipse (68% CL and 90% CL contours are shown) is consistent with the banana-shaped 
region of the S-T plane ahowed by the standard model, so long as the Higgs mass is small. The 
fit result is presently consistent with S and T = 0. Also shown are the itlcj constraints for four 
representative electroweak observables: sin^ O^y, M^, Tz, and R~ (from neutrino deep-inelastic 
scattering), which appear as bands with different slopes. At the present level of precision in the 
global electroweak data, the S-T analysis has excluded a subclass of the models that propose 
strong symmetry breaking mechanisms [conventional "technicolor" models (^)] that predict 
significant deviations from S = 0. The S-T analysis is beginning to probe hypothetical extensions 
of the standard model, such as super symmetry. For example, the region of the S-T plane allowed 
by the minimal supersymmetric model (MSSM) overlaps the light Higgs end of the standard- 
model zone but also includes larger deviations from T = 0. Some authors have already explored 
how the precision electroweak data constrain the free parameters of the MSSM (pgj). 
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We have reviewed highlights of the Z°-pole physics program carried out by the SLD collabo- 
ration at the SLC. The combination of beam polarization, small beam sizes, and a CCD pixel 
vertex detector has led to a series of measurements that complement those performed at Fer- 

milab and CERN and that have probed the standard model and its extensions. Many of these 
measurements provide unique insight into the physics of fundamental interactions. In addition, 
the novel techniques developed for these measurements have applicability to future colliding 
beam experiments. 

An important result of the precision electroweak program at SLD is that a standard-model 
Higgs boson, if it exists, should be relatively light. This conclusion is also supported by a fit 
to the complete set of worldwide precision electroweak data. In particular, these fits indicate 
that the Higgs boson may be within reach of the Tevatron experiments at Fermilab in the next 
several years. If the Higgs boson eludes the Tevatron experiments, a much greater reach is 
expected once the CERN LHC program begins in about 2006. Also important are the SLD 
electroweak measurements performed with individual lepton and quark flavors. This diverse 
set of measurements clearly illustrates the versatility of the powerful flavor-separation methods 
available to an experiment at a linear collider. Finally, studies of B^-B^^ mixing at SLD have 
contributed significantly to the worldwide effort to constrain the quark mixing matrix and CP 
violation in the standard model. These played an important role, complementary to the CP 
violation measurements currently under way. 

Looking a bit further ahead, high-energy (~500 GeV and up) linear collider projects are 
presently under consideration in the United States, Germany, and Japan. These new collider 
projects would build on the experience gained from the construction and operation of the SLC 
and SLD, as well as from the physics results and techniques developed there. 
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